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Abstract 

The problem of metastability for a stochastic dynamics with a parallel updating rule is ad- 
dressed in the Freidlin-Wentzel regime, namely, finite volume, small magnetic field, and small 
temperature. The model is characterized by the existence of many fixed points and cyclic pairs 
of the zero temperature dynamics, in which the system can be trapped in its way to the stable 
phase. Our strategy is based on recent powerful approaches, not needing a complete description 
of the fixed points of the dynamics, but relying on few model dependent results. We compute 
the exit time, in the sense of logarithmic equivalence, and characterize the critical droplet that 
is necessarily visited by the system during its excursion from the metastable to the stable state. 
We need to supply two model dependent inputs: (1) the communication energy, that is the 
minimal energy barrier that the system must overcome to reach the stable state starting from 
the metastable one; (2) a recurrence property stating that for any configuration different from 
the metastable state there exists a path, starting from such a configuration and reaching a lower 
energy state, such that its maximal energy is lower than the communication energy. 
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1. Introduction 

Metastable states are very common in nature and are typical of systems close to a first 
order phase transition. It is often observed that a system can persist for a long period 
of time in a phase which is not the one favored by the thermodynamic parameters; clas- 
sical examples are the super-saturated vapor and the magnetic hysteresis. The rigorous 
description of this phenomenon in the framework of well defined mathematical models 
is relatively recent, dating back to the pioneering paper [CGOV], and has experienced 
substantial progress in the last decade. See [OV] for a list of the most important papers 
on this subject. 

A natural setup in which the phenomenon of metastability can be studied is that 
of Markov chains, or Markov processes, describing the time evolution of a statistical 
mechanical system. Think for instance to a stochastic lattice spin system. In this context 
powerful theories (see [BEGM, MNOS, OS]) have been developed with the aim to find 
answers valid with maximal generality and to reduce to a minimum the number of model 
dependent inputs necessary to describe the metastable behavior of the system. Whatever 
approach is chosen, the key model dependent question is the computation of the minimal 
energy barrier, called communication energy, to be overcome by a path connecting the 
metastable to the stable state. Such a problem is in general quite complicated and becomes 
particularly difficult when the dynamics has a parallel character. Indeed, if simultaneous 
updates are allowed on the lattice, then no constraint on the structure of the trajectories 
in the configuration space is imposed. Therefore, to compute the communication energy, 
one must take into account all the possible transitions in the configuration space. 

The problem of the computation of the communication energy in a parallel dynamics 
setup has been addressed in [C, CN]. In particular, in [CN] the typical questions of 
metastability, that is the determination of the exit time and of the exit tube, have been 
answered for a reversible Probabilistic Cellular Automaton (see [GLD,R,St,To,D,CNP]), 
in which each spin is coupled only with its nearest neighbors. In that paper it has been 
shown that, during the transition from the metastable minus state to the stable plus 
state, the system visits an intermediate chessboard-like phase. In the present paper we 
study the reversible PCA in which each spin interacts both with itself and with its nearest 
neighbors; the metastable behavior of such a model has been investigated on heuristic and 
numerical grounds in [BCLS]. The addition of the self-interaction changes completely the 
metastability scenario; in particular we show that the chessboard-like phase plays no role 
in the exit from the metastable phase. 

Another very interesting feature of this model is the presence of a large number of fixed 
points of the zero-temperature dynamics in which the system can be trapped. Following 
the powerful approach of [MNOS], we can compute the exit time avoiding a complete 
description of the trapping states. However, we cannot describe the exit tube, i.e., the 
tube of trajectories followed by the system during its exit from the metastable to the 
stable phase. The only information on the exit path that we prove in this paper is the 
existence of a particular set of configurations which is necessarily visited by the system 
during its excursion from the metastable to the stable state. This set plays the role of 
the saddle configuration set, which is usually introduced in the study of the metastable 
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behavior of sequential dynamics. 

According to the approach of [MNOS], the model dependent ingredients that must 
be provided are essentially two: (1) the solution of the global variational problem for all 
the paths connecting the metastable and the stable state, i.e., the computation of the 
communication energy; (2) a sort of recurrence property stating that, starting from each 
configuration different from the metastable and the stable state, it is possible to reach a 
configuration at lower energy following a path with an energy cost strictly smaller than 
the communication energy. 

To solve the global variational problem (see items [2] and [3] in Theorem 12. 3ft . we ob- 
tain an upper bound on the communication energy by exhibiting a path connecting the 
metastable state to the stable state whose maximal energy is equal to the communication 
energy. To find the lower bound, we perform a partition of the configuration space, study 
the transitions between configurations in these partitions, and reduce the computation to 
the optimal one (see Figured]). To prove the recurrence property (see item [T] in Theo- 
rem we have to face the problem of the existence of a large number of fixed points 
of the dynamics. We solve this problem by showing that, for each configuration different 
from the metastable state, it is possible to find a path connecting it to the stable state, 

1. e., to the unique global minimum of the energy, such that the energy along this path is 
strictly smaller than the communication energy. 

We finally give a brief description of the content of the paper. In Section [2] we define 
the model and state our main result in Theorem 12.11 The proof of Theorem I2.1[ based 
on the model dependent results in Theorem 12.31 and on [MNOS], is given in Section [2781 
Section [3] is devoted to the proof of the estimates on the energy landscape stated in 
Theorem 12.31 namely, the global variational problem (items [2] and [3]) and the recurrence 
property (item[T]). The proof of items [2] and [3] relies on Proposition 13.21 which is proven 
in Section [H The appendix is devoted to a brief review of results in [MNOS]. 

2. Model and results 

In this section we introduce the basic notation, define the model, and state our main result. 
In particular, Sections I2.1H2.4I are devoted to the definition of the Probabilistic Cellular 
Automaton which will be studied in the sequel. In Section 12.51 we state Theorem 12.11 
with the results on the metastable behavior of the system. In Section 12.61 we introduce 
the transition rates and the zero temperature dynamics; in Section 12.71 we develop an 
heuristic argument on which the proof of the theorem is based. In Section 12.81 finally, we 
prove Theorem 12.11 

2.1. The lattice 

The spatial structure is modeled by the two-dimensional finite square A := {0, . . . , L — l} 2 , 
where L is a positive integer, with periodic boundary conditions; note that A is a torus. 
We shall use the metric induced by the Euclidean distance on the flat torus. An element 
of A is called a site. We use X c := A \ X to denote the complement of X C A. 

Let x G A; we say that y G A is a nearest neighbor of x if and only if the distance on 
the torus of x from y is equal to 1. For X C A, we say that y G X c is an element of the 
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external boundary dX of X if and only if at least one of its nearest neighbors belongs to 
X; we let also X := X U dX be the closure of X. Two sets X, Y C A are said to be not 
interacting if and only if for any i6l and their distance on the torus is larger or 

equal to a/5. 

Let x = (xi,X2) G A; for £1,^2 positive integers we let Qe lt £ 2 (x) be the collection of 
the sites ((a;i + Si) mod L, (x 2 + s 2 ) mod L) for Sj = 0, . . . , Xi + l{ — 1 where i = 1,2. 
Roughly speaking, Qe lt e 2 {x) is the rectangle on the torus of side lengths t\ and £ 2 drawn 
starting from x and moving in the positive direction along the two coordinate axes. For 
£ a positive integer we let Qe(x) := Qe/(x). 

2.2. The configuration space 

The single spin state space is given by the finite set {—1, +1}; the configuration space in 
X C A is defined as Sx '■= { — 1, +1} X and considered equipped with the discrete topology 
and the corresponding Borel a algebra Tx- The model and the related quantities that 
will be introduced later on will all depend on A, but since A is fixed it will be dropped 
from the notation; in this spirit we let <Sa =: S and T\ =: T . 

Given a configuration a E S and X C A, we denote by ox the restriction of a to X. 
Let m be a positive integer and let X\, . . . , X m C A be pairwise disjoint subsets of A; for 
o"fe G Sx k , with k = 1, . . . , m, we denote by <7ier 2 • ■ ■ a m the configuration in Sxiu-ux m 
such that (o~icr2 ■ • • o~m)x k — o~k for all k G {1, ... , m}. Moreover, given a G S and iG A, 
we denote by a x the configuration such that a x (x) = — o~(x) and cr x (y) = o~(y) for y ^ x. 
Let x G A, we define the shift Q x acting on S by setting (Q x o~) y := a y+x for all y G A and 
a eS. 

Given a function / : S — > M, if / G J^x we shall sometimes write f(&x) f° r /( (7 )- Let 
f,g:S^Sbe two functions, we consider the product or composed function fg:S^S 
such that fg(cr) := f(g(cr)) for any o~ E S. We also let / 2 := // and, for n a positive 
integer, f n := ff n . We say that a configuration a G 5 is a fixed point for the map 
/ : S — > S if and only if /(er) = er. Let cr G 5, consider the sequence f n (o') with n > 1, if 
there exists n' such that f n {c) = f n (a) for any n > n', we then let /cr := /" cr. 

2.3. The model 

Let /5 > and h E R such that < 1 and 2/ his not integer. We consider the Markov 
chain on 5 with transition matrix 



P(<7>V) '■= (v( x )) Vo-,?7G5 (2.1) 



where, for each x E A and a E S, p x ,a{') is the probability measure on S{ x y defined as 
follows 

P*,«r(s) := t- 7 9/ o\ Q 7 r—TTT = \ [1 + s tanh (3 (S a (x) + h)\ (2.2) 

1 + exp {— 2ps{b a [x) + h)\ 2 

with s G { — 1, +1} and 

:= Y^iy) (2.3) 
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The normalization condition p Xj(T (s)+p X)(T (—s) = 1 is trivially satisfied. Note that p x ,-(s) £ 
J~rpr for any x and s, that is the probability p x ,a( s ) f° r the spin at site x to be equal to s 

depends only on the values of the five spins of a inside the cross {x} centered at x. 

Such a Markov chain on the finite space S is an example of reversible probabilistic 
cellular automata (PCA), see [D,CNP]. Let n £ N be the discrete time variable and let 
o~ n £ S denote the state of the chain at time n; the configuration at time n + 1 is chosen 
according to the law p(a n , •), see (12.11) . hence all the spins are updated simultaneously and 
independently at any time. Finally, given a £ S we consider the chain with initial config- 
uration <r = a, we denote with F a the probability measure on the space of trajectories, 
by IE CT the corresponding expectation value, and by 

r a A := inf{t > : a t £ A} (2.4) 

the first hitting time on A C S. We shall drop the initial configuration from the notation 
(I2.4jl whenever it is equal to — 1, i.e., we shall write ta instead of t^~. 

2.4- The stationary measure and the phase diagram 

The model ( 12. II) has been studied numerically in [BCLS]; we refer to that paper for 
a detailed discussion about its stationary properties. Here we simply recall the main 
features. It is straightforward, see for instance [CNP,D], that the PCA (12. ip is reversible 
with respect to the finite volume Gibbs measure /i(e>") := exp{— H(a)}/Z with Z : = 

E^ e s ex P{-- fr ( 7 ?)} and 

H(a) := H p>h (a) := -/3h £ a(x) - £ log cosh \J3 (S a (x) + h)} (2.5) 

In other words the detailed balance condition 

p(a, 7]) e- H(a) = p{r], a) e~ H{r]) (2.6) 

is satisfied for any a, rj £ S; hence, the measure /x is stationary for the PCA ( 12. II) . In order 
to understand its most important features, it is useful to study the related Hamiltonian. 
Since the Hamiltonian has the form ( 12. 5p . we shall often refer to 1/(3 as to the temperature 
and to h as to the magnetic field. 

The interaction is short range and it is possible to extract the potentials; following 
[BCLS] we rewrite the Hamiltonian as 

H p>h {a) = ^ U xAh {a) -PhJ2 <x) ( 2 -7) 

2 e A xeA 

where U x ^^{o') = Uq^^{^>xO') 1 recall that the shift operator Q x has been defined in 
Section |2~21 and that periodic boundary are considered on A, and 

uo Ah {°) = - Yi j w> h n a ^ ( 2 - 8 ) 

xcjo} x&x 
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The six coefficients Jo,p,h, ■ ■ ■ > Jb,p,h are determined by using (12. 5p . (12. 7p . and (12. 8p . In 
the case h = only even values of |X| occur and we find that the pair interactions are 
ferromagnetic while the four-spin interactions are not. For a more detailed discussion 
see [BCLS]. 

The definition of ground state is not completely trivial in our model, indeed the Hamil- 
tonian H depends on (3. The ground states are those configurations on which the Gibbs 
measure /i is concentrated when the limit (3 — > oo is considered, so that they can be 
defined as the minima of the energy 



Let X C iS, if the energy E is constant on X, we shall misuse the notation by denoting 
by E(X) the energy of the configurations in X . 

We first consider the case h — 0. Since E(a) = — XLeA \Sct{x)\, h is obvious that 
there exist the two minima +1, —1 G S, with ±I(x) = ±1 for each x G A, such that 
E(+l) = E(-l) = -5|A|. For h ^ we have E(+l) = -\A\{h + |5 + h\) and E(-l) = 
—\A\(—h + | - 5 + h\); it is immediate to verify that E(+l) < E(—l) for h > and 
E(—l) < E(+l) for h < 0. We conclude that at h — there exist the two ground states 
— 1 and +L At /i > the unique ground state is given by +1 and at h < the unique 
ground state is given by —1. The phase diagram at finite large (3 and h = has been 
studied rigorously in [DLR]. 

2.5. Metastable behavior 

We pose now the problem of metastability and state the related theorem on the exit time. 
In this context, configurations with all the spins equal to minus one excepted those in 
rectangular subsets of the lattice will play a key role. We then let 



for any a G S; the set A + (a) will be called the support of a. We say that a G S is a 
rectangular droplet with side lengths £ and m, with £, m integers such that 2 < £, m < L—2, 
if and only if there exists x G A such that either A + (o") = Qe >m (x) or A + (o") = Q m ,e(x). 
We say that a G 5 is a n-rectangular droplet with side lengths £x,mx, . . . ,£ n ,m n , with 
n > 1 an integer and £i, rrii integers such that 2 < £i, rrii < L — 2 for z = 1, . . . , n, if and 
only if A + (<t) is the union of n pairwise not interacting rectangles (see Section I2TT1) with 
side lengths £i and mi for i — 1, . . . , n. We finally say that a G S is a multi-rectangular 
droplet if and only if a is a n-rectangular droplet for some integer n > 1. Note that a 1- 
rectangular droplet is indeed a rectangular droplet. Square droplets are defined similarly. 

Consider, now, the model (12. ip with < h < 1 and suppose that the system is 
prepared in the state o"o = —1; in the infinite time limit the system tends to the phase 
with positive magnetization. We shall show that the minus one state is metastable in 
the sense that the system spends a huge amount of time close to —1 before visiting +1; 
more precisely the first hitting time r+i to +1 (recall (I2.4p and the remark below) is an 
exponential random variable with mean exponentially large in (3. 
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(2.9) 



(a) := {xGA: a{x) = ±1} 



(2.10) 



Moreover, we give some information on the exit path that the system follows during 
the escape from minus one to plus one. More precisely, we show that there exists a class of 
configurations C G S, called set of critical droplets, which is visited with high probability 
by the system during its escape from — 1 to +1. Let the critical length A be defined as 



A := 



(2.11) 



where, for any positive real x, we denote by [^J the integer part of x, i.e., the largest 
integer smaller than or equal to x. Since h has been chosen such that 2/ his not integer, 
see Section 12.31 we have that A = 2/h + 5h with 5h G (0, 1). The set C is defined as the 
collection of configurations with all the spins equal to —1 excepted those in a rectangle 
of sides A — 1 and A and in a pair of neighboring sites adjacent to one of the longer sides 
of the rectangle. 

Given 7 G C we let 

T := E(i) - E(-l) + 2(1 + h) = -4h\ 2 + 16A + 4h(X - 2) + 2(1 + h) (2.12) 

Note that by fl2~T2l and (12TTTT) it follows 

T<8A + 10-2/i (2.13) 

The simple bound above will be used in Section 13.21 to prove (I3.25P and in Section 14.41 

As has been explained in the introduction, the energy of the configurations in C is 
strictly connected to the typical exit time from the metastable state, indeed we have the 
following theorem. 

Theorem 2.1 For h > small enough and L = L(h) large enough, we have that 
1. for any e > 

lim P_ 1 (e /?r -^ £ < r+i < e^ +/3e ) = 1 (2.14) 

2. 

lim ilogE_ 1 [r +1 ]=r (2.15) 

p— >oo p 



3. 



lim P_i(r c < 7+1) = 1 (2.16) 

P~^oo 



In other words, the above theorem states that the random variable (1//3) log 7+1 converges 
in probability to V as (3 — ► 00 and that the logarithm of the mean value of r + i divided 
times j3 converges to V in the same limit. Moreover, the last item ensures that, before 
reaching the stable state +1, the system started at — 1 must necessarily visit the set of 
critical droplets C. 

The proof of Theorem 12.11 will be given in Section 12.81 We note that, as usual in 
Probabilistic Cellular Automata (see also [CN]), the highest energy V reached along the 
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e -2/3(5+h) e - 2 /3(5-h) 



e -2/3(3+h) e - 2 /3(3-h) 

+ 

2 f3(l+h) 2 p(l-h) 



+ 



1 _ £-2/3(1-/*) - - + 1 - e~ 2/3(1+h) 

+ 

1 _ p -2/3(3-h) + 1 _ p -2f3(3+h) 

+ 

Figure 1: Large j3 behavior of the probabilities for the flip of the central spin for all 
possible configurations in the 5-spin neighborhood. 

exit path is not achieved in a configuration, which is the typical situation in Glauber 
dynamics. Such a V is the transition energy (see definition (I2.18P ) of the jump from 
the "largest subcritical" configuration to the "smallest supercritical" one, see also the 
heuristic discussion in Section |2~71 

2.6. Transition rate and zero temperature dynamics 

In our problem (see also [CN]) the energy difference between two configurations a and 77 
is not sufficient to establish whether the system prefers to jump from o to 77 or vice versa. 
Indeed, for some pairs of configurations a sort of barrier is seen in both directions; more 
precisely, it is possible to find a and 77 such that both p(o~, 77) and p(rj, a) tend to zero 
in the zero temperature limit (3 — > 00. As an example of such a behavior, consider the 
two following configurations: a is such that all the spins are equal to minus one excepted 
those associated with the sites belonging to an I x I rectangle, with 3 < i < L — 2, and to 
a two-site protuberance attached to one of the sides of the rectangle; 77 is a configuration 
obtained starting from a and flipping the spin associated with one of the sites neighboring 
both the rectangle and the protuberance. By using fl2.1|) - fl2.3p . it is easy to show that 
p(<7, 77) ~ exp{— 2/5(1 — h)} and p(r], a) ~ exp{— 2/5(1 + h)} for large /3; see also Figured! 
where we have reproduced the table in [BCLS, FIG. 1] with the list of the single site event 
probabilities. In that figure, the large (3 behavior of the probability, associated to the flip 
of the spin at the center, is computed. 

To manage those barriers we associate the transition Hamiltonian H(a, 77) to each pair 
of configurations cr, n e S. More precisely we extend the Hamiltonian (12.51) to the function 
H:SUSxS^Rso that 

H(a,f)) :=H{a)-\ogp{a,n) (2.17) 
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By the detailed balance principle (12. 6p . we have H(a, rf) = H(r), a) for any a,rj G S. Note 
that, by definition, if (17,77) > max{if(er), H{rf)} and 77) = exp{ — [H(a, rf) — H(a)}}; it 
is then reasonable to think to H(a, rf) as to the Hamiltonian level reached in the transition 
from a to 77. As already noted in Section 12.41 since the Hamiltonian depends on j3, it is 
useful to compute its limiting behavior. We then define the transition energy 

E(a,7j):= lim ^H(a, V ) (2.18) 

Note that by using f]2.9f) . ( 12.171) . and the symmetry of the transition hamiltonian, we get 
E(o~, rf) = E(j], a) and E(a, rf) = E(a) + A(a, rf) > max{E(a), E{rf)} (2.19) 
for any o~,rj E S, with 

A(tr, rf) := - lim - \ogp(a, rf) = V 2\S (x) + h\>0 (2.20) 

Ti(x)(So-{x)+h)<Q 

the transition rate] notice that in the second equality we have used the definition (12. 9p of 
E(a), (EI7|), jZU), and (Q. 

The non-negative transition rate A will play a crucial role in the study of the low 
temperature dynamics of the model (12.11) ; indeed it can be proven that the model satisfies 
the FW condition in [OV, Chapter 6], that is for any a, r] e S and j3 > large enough 

e -pA(*, v )-f3 7 (p) < p ^ ^ < e -m^n)+PiW (2.21) 

where 7(/3) does not depend on a, 77 and tends to zero in the limit (3 — > 00. From (12.211) 
it follows that p(o~,r)) — > 1 for /3 — * 00 if and only if A(cr, 77) = 0. On the other hand, 
if A(ct, 77) > 0, then p(a,rj) — > exponentially fast and with rate A(cr, 77) in the limit 
/3 — > 00, so that A can be interpreted as the cost of the transition from a to 77. 
To get (12.211) we first prove that for (3 large enough 

-^[H(a,rj)-H(a)} + [E(a,ri)-E(a)}\<e-^ (2.22) 

The bound fl2T2T|) shall follow easily from fl2T22|) . (j2TT7j) . and the second equality in fl2TT9l) 
relating the transition energy to the transition rate. To prove (I2.22p we note that by using 
(J2TZD, (J2U), O, and dOHD we get 

-\H(a, 77) - //(a)] - [E(a, 77) - = ± £ log(l + e -^W+*l) (2.23) 
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^[H(o,r))-H(o)]-[E(o,r])-E(o)} 



i Ml + e- 2 ^ x)[s ° {x)+h] ) + ^ 2 V (x)[S a (x) + fc] 

77(a;)(S CT (a;)+fe)<0 

- ^ log(l + e - 2 Pv(x)[SAx)+h]^ + I log ^ + e -2/3r ? (x)[5 CT (x)+h]^ 

^ 16A: ^ ISA: 

77(x)(S CT (a;)+/i)>0 77(i)(S CT (x) + h)<0 

+ ^ 2 V (x)[S a (x) + h] 



r,(i)(S CT (i) + h)<0 

= 1 l og (l + e -m(*)[S„(x)+h]j + I log( e +2j8»j(*)[&(*)+fc] + J) 

^ i6A: ^ i6A: 

r;(x)CS CT (x)+h)>0 >i(i)(S„(i)+)i)<0 

yielding (EZSj) . The bound 112T22D follows once we note that log(l+exp{-2/5|5 a (x)+/i|}) > 
for any x G A and |5 CT (a;) + /i| > 1 — /i uniformly in <r G 5 and x G A, and choose 
0>(log|A|)/(l-/i). 

We finally introduce the zero temperature dynamics. Consider a configuration cr G S 
and s G {— 1,+1}; since < 1, from 1 12. 2ft it follows that the probability p x ,a(s) tends 
either to or to 1 in the limit (3 — > oo. Thus, due to the product structure of 1 12.11) . given 
cr there exists a unique configuration r] such that p(o,rj) — > 1 in the limit /3 — > oo. This 
configuration is the one such that each spin rj(x) is chosen so that p x .a{v( x )) ~^ 1 f° r 
j3 — > oo. We introduce the map T : S — »• S, called the zero temperature dynamics, which 
associates to each o G S the unique configuration To such that p(o, To) — ► 1 in the limit 
j3 — > oo. 

Lemma 2.2 Given o,r] E S, we have that A(o,rj) = if and only ifrj — To. 

Proof of Lemma \2.2\ The lemma follows immediately by using the definition of the zero 
temperature dynamics T and the remarks below 112. 21[) . □ 



2. 7. Stable states and stable pairs 

The proof of Theorem 12 .11 although mathematically complicated, relies on a very straight- 
forward physical argument based on a careful description of the low temperature, i.e., large 
(3, dynamics. In this section we give an heuristic explanation of the exponential estimate 
112. 14p for the exit time r+i. 

We introduce, first, the notion of stable configurations. If To = o the configuration o 
is called stable; equivalently, we say that o G S is stable if and only if for any r\ G S \ {o} 
one has p(o, rj) — > in the limit (3 — > oo. If o is not stable and T 2 o = o, we say that 
(cr, To) is the stable pair associated to cr, equivalently we say that (cr, TV) is a stable pair if 
and only if p(o, To) — > 1 and p{To, o) — > 1 in the limit /3 — ■> oo. Recall A is non-negative, 
by (12. 1 9j) and Lemma [2.21 it follows 

E(o,To) = E(o) and E(o) > E(To) (2.24) 
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Figure 2: Examples of stable states, pluses and minuses are represented respectively by 
grey and white regions. 



for any a G S. 

Note that a stable pair (0, 77) is a 2-cycle of the map T, indeed Ta = r\ and T?7 = a. It 
is easy to show that cycles longer than two do not exist for such a map. Suppose, by the 
way of contradiction, that cr 1; . . . , o~ n G S, with n > 3 integer, are such that 0, 7^ <7j for 
i ^ j, T<Ji = a i+ i for % — 1, . . . , n — 1, and Ta„ = ui. By the inequality in (12.241) . it follows 
that E(a\) > ••• > E(a n ) > i?(0"i), which implies -E(0i) = ••■ = E(a n ). This result, 
together with the equality in (12.241) and (12.191) . implies that A(o"i,o" 2 ) = A(cr 1 ,cr n ) = 0. 
Hence, by recalling Lemma T2.2} we get To~\ = 02 and Tax — &n- By definition of the map 
T, we finally get o n = <J2, which contradicts the hypothesis that <Ji ^ <jj for i 7^ j. 

As mentioned above, our model is characterized by the presence of a large number of 
stable configurations. Indeed, only those configurations in which there exists at least one 
spin with a majority of opposite spins among its neighbors are not stable, see Figured] All 
the configurations in which each spin is surrounded by at least two spins of the same sign 
are, instead, stable; some of the possible situations are shown in Figure [2 In particular, 
notice that plus squared rings plunged into the sea of minuses are stable states. This 
scenario is complicated by the presence of stable pairs; some of them are depicted in 
Figure [3j Notice, in particular, the chessboards leaned to stable pluses regions. As we 
shall see in the sequel, the stable pairs do not play any important role in the study of 
metastability in model (12. ip . We also recall that, in the case of a similar model studied 
in [CN], due to the presence of such pairs, the system was forced to visit an intermediate 
chessboard phase in its way from the minus metastable phase to the stable plus phase. 

We describe, now, the typical low temperature behavior of the dynamics. Suppose 
that the initial condition is a = a G S; at low temperature, with high probability, the 
system follows the unique zero temperature trajectory 

0- o = a, 0i = Ta, 2 = Ta x = T 2 a, . . . , a t = T(T*~V) = TV, . . . 
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Figure 3: Examples of stable pairs, pluses and minuses are represented respectively by 
grey and white regions. 



Once the zero temperature trajectory ends up in a stable configuration, it remains there 
forever. Different trajectories are observed with probability exponentially small in (3. 

We can now depict the typical behavior of the system at very low temperature. Recall 
the definitions given in the last paragraph of Section [2721 Starting from a, the system will 
reach in a time of order one either the stable configuration Ta or the stable pair associated 
to T 2 a; note that Ta and T 2 a are unique. After a time exponentially large in (3, the chain 
will depart from the stable configuration, or from the stable pair, and possibly reach a 
different stable configuration, where it will remain for another exponentially long time. 
And so on. It is then clear that, in the study of the low temperature dynamics, a key role 
is played by stable configurations and stable pairs; indeed a large amount of the time of 
each trajectory is spent there. 

Among the large number of possible stable states, there are those configurations in 
which the plus spins fill a rectangular region; recall the definition of rectangular droplets 
given at the beginning of Section 12.51 In [BCLS] it has been conjectured that those 
rectangular stable configurations are the relevant ones for metastability. Moreover, there 
has been developed an heuristic argument to show that A, see (12. lip , is the critical length 
in the sense explained below. Rectangular droplets with smallest side length smaller or 
equal to A — 1 are subcritical, namely, starting from such a configuration the system visits 
— 1 before +1 with probability tending to one in the limit (3 — > oo. Rectangular droplets 
with smallest side length larger or equal to A are supercritical, namely, starting from such 
a configuration the system visits +1 before —1 with probability tending to one in the 
limit j3 — > oo. 

We reproduce shortly the heuristic argument in [BCLS, Section IV] yielding the above 
conclusions. Consider a square droplet of side length £; we shall identify the best growth 
and shrinking mechanisms and, by comparing the related typical times, get the critical 
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Figure 4: Shrinking mechanism. 



length. First note that the configuration obtained by attaching a single site protuberance 
to one of the sides of the droplet is not stable (see Figured]); it is needed at least a two-site 
protuberance to get a stable configuration. The parallel dynamics allows the formation of 
a two-site protuberance in one step; from Figure [Hand the product structure of ( 12. ip . it 
follows that the typical time for this process is r onc ~ exp{4/3(3 — h)}. On the other hand, 
the protuberance can be formed in two consecutive steps: first a single site protuberance 
appears and, then, one of the two minuses adjacent both to the square droplet and to the 
protuberance is flipped. By using again the data in Figure [I] we get that the typical time 
for the two-step process is r two ~ exp{2/3(3 — h) + 4/3(1 — h)}, where 2/3(3 — h) is the 
cost of the first step and 4/3(1 — h) is the sum of the costs paid in the second step to keep 
the single site protuberance and to flip the adjacent spin. Clearly r two <C r onc for (3 large; 
hence, the most efficient mechanism to produce a two-site protuberance is the two-step 
one. 

The presence of a two-site protuberance is sufficient to ensure the growth of the 
droplet. Indeed, noted that exp{2/3(l — h)} is the smallest typical time needed to leave a 
stable configuration (see Figured]), we have that the side with the two-site protuberance 
is filled by pluses via a sequence of £ — 2 flips of a minus spin with two neighboring pluses. 
Since each of those flips happens with typical time of order exp{2/3(l — h)}, we conclude 
that the growth time Tgr OW th is equal to r two . 

For what concerns the shrinking mechanism, it is easy to show that the most efficient 
one is the flipping of plus spins having two neighboring minuses (corner erosion). The 
shrinking is then performed via a sequence of configurations as in Figure H] requiring the 
erosion of i— 1 corner pluses and the final flipping of the unstable single site protuberance. 
Note that the intermediate configurations, joining the starting I x i square droplet to the 
ending single protuberance configuration, are stable; their lifetime, i.e., the typical time 
that must be waited for to see the system performing a transition, is exp{2/3(l — h)} 
(see Figured]). It follows that suitably long persistence in the 1 — 1 intermediate stable 
configurations must be provided for in the most efficient shrinking path. The rate at which 
the entire process occurs is thus estimated as the rate for one erosion, exp{— 2/3(1 + h)}, 
times the probability that I — 2 further erosions occur within the lifetime exp{2/3(l — h)}, 
which is of order [exp{— 2j3(l+h)} exp{2/3(l— h)}] e ~ 2 . We then conclude that the shrinking 
time is estimated as r shrin king ~ exp{2^(l + h) + (£ - 2) [2/3(1 + h) - 2/3(1 - h)]}. 

By comparing, finally, r s h r inking and r growt h, we get that growth is favored w.r.t. shrink- 
ing if and only if I > [2/h\ +1. This remark strongly suggests that the length A, defined 
in (12. lip , plays the role of the critical length for what concerns the metastable behavior 
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of the model. 

We come, finally, to the heuristic argument suggesting the estimate (12.141) for the 
exit time. It is reasonable to suppose that the exit path visits an increasing sequence of 
subcritical rectangular droplets, whose side lengths differ at most by one. The highest 
energy along such a path wil be attained in the segment leading from the largest subcritical 
A x (A — 1) droplet to the smallest supercritical A x A droplet. More precisely, denote 
by 7r the configuration obtained by attaching a single site protuberance to one of the two 
longer sides of the A x (A — 1) droplet and by 7 the configuration obtained by nipping 
in 7r a minus spin adjacent to the rectangle and neighboring the single site protuberance. 
Recall the discussion above about the growth mechanism. It follows that the highest 
energy along the exit path must be attained in the transition from it to 7, so that it is 
equal to E(ir,j) (see (I2.18P ). It is then reasonable to expect that the typical exit time is 
of order exp{P[E(n, 7) — E(— 1)]}. Using the expression 

E(i/>) - E(-l) = -4M x l 2 + 8(^ + £ 2 ) (2.25) 

for a rectangular droplet ip G S of side lengths £\ and £2, recall that in such a configuration 
2 < £\, £ 2 < L — 2, it is an easy exercise to show that E(ir, 7) — E(—V) = T, see (12.121) . 

2.8. Escape time 

In this section we prove Theorem 12.11 The main ingredients will be the general results 
[MNOS, Theorem 4.1, 4.9, and 5.4], the solution of the model dependent variational 
problem (12.291) . i.e., the computation of the energy barrier between — 1 and +1, and the 
recurrence estimate (12.281) . In [MNOS] the theory has been developed with quite strict 
hypotheses on the dynamics, see [MNOS, equation (1.3)], nevertheless it can be shown 
that the same results hold in the present setup, see Appendix [A] 

To state the estimates on the energy landscape we need few more definitions. A finite 
sequence of configurations uj = {oui, . . . ,u> n } is called path with starting configuration uj 1 
and ending configuration uj n ; we let |u;| := n. We let ft, := (S N \{°} be the collection of 
all the possible paths. Given two paths uj and u', such that u\ u \ = uj[, we let u + uj' : = 
{loi, . . . , u)\ u \,L)2, . . . , wL/i}; note that \uj + u/| = |a>| + \oj'\ — 1. Given a path to, we define 
the height along uj as 

$ := / EM if M = 1 (2.26) 

[ maXj = i r .. i | a; |_ 1 E(uji, u^+i) otherwise 

Let A, A' C S, we denote by Q(A,A') the set of all the paths uj G Q such that uji G A 
and u\ u \ G A', that is the set of paths starting from a configuration in A and ending in a 
configuration in A'. The communication energy between A, A' C S is defined as 

$(A,A'):= min $ w (2.27) 

lo£0(A,A') 

If A = {cr}, we shall misuse the notation by writing 0(cr, A') instead of Q({a},A') and 
$(<7, A') instead of $({cr}, A'). 
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Theorem 2.3 Recall the definition of T in A2.12\) . Suppose that h > is chosen small 
enough. Then 



1. for any a G S \ { — 1} 



$(a, +1) - E(a) < T 



(2.28) 



H-l, +i)-E(-l) = T 



(2.29) 



3. for each path uj = {oui, . . . , u n } G 0(— 1, +1) such that $ w — E{— 1) = V , there exists 
i G {2, . . . , n} such that Ui G C and E(u)i-i,Ui) — E(—l) = T. 

Theorem 12.31 will be proved in Section [3J Recall Theorems IA.lHA.3l and the definitions 
given before them. 

Proof of Theorem \2.1[ By using the results discussed at the end of Section 12.41 we 
have S s = {+!}. We remark that, since 5 s = {+!}, then for any a G S \ S s we have 
E(+l) < E(a); this implies, together with (12^91) and ([ZZ8D , that S m = {-1} and 
V-i = T. Finally, items [U and [2] follow from Theorems IA.1I and IA.21 respectively. 

Proof of item [3J By using item [3] in Theorem 12.31 we get that C is a gate for the 
transition from — 1_ to +1_. The item [3] in Theorem 12.11 follows from Theorem IA.3I □ 



3. The recurrence property and the variational problem 

In this section we prove the energy landscape estimates stated in Theorem 12. 3( in partic- 
ular the recurrence property fj2.28j) is proven in Section 13.11 and the variational problem 
f)2.29p is solved in Section 13.21 The proof of items [2] and [3] of Theorem 12.31 relies on 
Proposition 13.21 which is stated in Section [3T21 and proven in Section HI 

We give few more definitions. Let a G S and x G A, we say that the site x is 
stable (resp. unstable) w.r.t. a if and only if a(x)S a (x) > (resp. a(x)S a (x) < 0). Note 
that the stable sites are those that are not changed by the zero temperature dynamics, 
more precisely Ta(x) = cr(x) if and only if x is stable w.r.t. a. Given a G S and 
k G {—5, —3, —1, +1, +3, +5} we denote by A fe (<r) the collection of the sites x G A ± (a) 
such that S a (x) = k, i.e., 



A±»:=A± 5 (<r)U-..UA±(<r) and A± fc (<r) := Af(tr) U ■ ■ • U A±(<r) (3.2) 



Finally, given a G S, we denote by A+(cr) (resp. A+(a)) the collection of the sites x G A 
such that cr(x) = +1 and x is stable (resp. unstable) w..r.t. a\ similarly we define A~(cr) 
and A~(<t). By definition of stable and unstable sites we get that, for any a G S, 



A+(a) = A+^a), A~(a) = A; +1 (a), A s +(a) = A+ +1 (a), and A s » = A^) (3.3) 



Af (a) := {X6A ± 



(a): S a (x) = k} 



(3.1) 



note that Ai 5 (er) = and A +5 (cr) = 0; moreover, we set 
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3.1. The recurrence property 

Equation (I2.28P in Theorem 12.31 states that, for any configuration a different from the 
metastable state —1, it is possible to exhibit a path uj joining a to the stable state +1, 
i.e., to the absolute minimum of the energy, such that ^ < E(o~) + T. On the heuristic 
ground, given a G S \ {—I}, there exists at least a plus spin; starting from such a plus it 
is possible to build a supercritical A x A droplet of pluses paying an energy cost strictly 
smaller than E(a) + Y. Indeed, by virtue of f)2.29p . starting from —1, the cost would be 
exactly V. On the other hand, starting from a, no energy must be paid to get the first 
plus spin and the other pluses of a, if any, help the production of the supercritical droplet. 

A rigorous proof needs the explicit construction of the path; such a path will firstly 
realize the growth of a supercritical A x A square with a as a background and then its 
growth towards +1. More precisely, recall A is a squared torus, let L be its side length 
and = (0, 0) the origin; recall the zero temperature dynamics mapping T defined in 
Section [2761 and let a G S be such that a(x) = +1 for any x G Q2,2(0). We define the path 



where the paths S n , with n = 2, . . . , L, and \l/ n , with n = 3, . . . , L, are constructed 
algorithmically. 

We first describe informally the algorithms. The path S n starts from the configuration 
£ n and ends in the configuration ip n+1 . The configuration £ n is such that the square 
Qn,n(0) = {0, . . . , n — 1} X {0, . . . , n — 1} is filled with pluses; the path S n fills with pluses 
the slice Qi, n (n, 0) = {(n, 0), . . . , (n, n — 1)}, adjacent to the square Qn !n (0), and produce 
ip n+1 in which the rectangle Q n +i, n (0) is filled with pluses. Similarly, the path starts 
from the configuration ip n and ends in the configuration £ n . The configuration ijj n is such 
that the rectangle Q n ,n-i(0) = {0, . . . , n — 1} x {0, . . . , n — 2} is filled with pluses; the path 
\l/ n fills with pluses the slice Q n ,i(0, n — 1) = {(0, n — 1), . . . , (n — 1, n — 1)}, adjacent to 
the rectangle Q n ,n-i(0), and produce £ n in which the square Q n ,n(0) is filled with pluses. 



Definition of S n . Let £ 2 := a, let n G {2, . . . , L — 1}, and suppose £ n is such that 
C(x) = +1 for x G <2n, re (0), then 

1. set z = 1, = T; 

2. if T 2 Q = Q then goto H else set i = i + 1 and Q = T^ n „ 1 and goto [21 

3. if £f(» = +1 for all x G Q 1>n (n, 0) then set tp n+1 = g and goto 

4. if Qi, n (n,0) n A+(£f) + 0, then pick y,y' G Qi,„(n,0) such that d(y,y') = 1, 

= -1, and y' G A+(ff), set i = i + 1, ^(y) = +1, Q^) = T&ifr) 
Wx G A \ {y} and goto |3j 

5. if Qi,„(n,0) n A+(£f) ^ 0, then pick y, y' G Qi,„(n,0) such that d(y,y') = 1, 

= -1, and ?/ G A+(tf), set z = i + 1, ^(y) = +1, OV) = +1, C(^) = 
T^ix) for any x G A \ {j/, y'} and goto EJ 
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n=3 



6. set % = i + 1, y = (n, 0), £f (y) = +1, (a) = T^^x) for any x G A \ {y} and 
goto El 

7. set ^ = i, E n = {£, ■ ■ ■ , and exit. 

At step [2] the algorithm follows the zero temperature dynamics down to the stable 
pair or to the stable state associated to £ n . At step [3] the algorithm checks if the slice 
Qi,n{ n :fy adjacent to the square Q n ,n(Q) is filled with pluses. In case of positive answer 
the algorithm jumps to step [7] and exits. If the answer is negative, then the slice is filled 
with pluses at steps |3H6] as follows: a minus adjacent to a stable plus is flipped (stepHJ); 
in absence of stable pluses, a minus adjacent to an unstable plus is flipped (step [5]). If the 
slice is filled with minuses, then the spin associated to the site (n, 0) is flipped (step [6]). 

Definition of Let n G {3, . . . , L} and suppose ip n is such that ip n {x) = +1 for 
x G Q„,n-i(0), then 

1. set i = 1, ip? = tfj n ; 

2. if T 2 ^ n = V>r then goto [3] else set i = i + 1 and = T^-i and goto El 

3. if ij%{x) = +1 for all x G Q n ,i(0,n - 1) then set £ n = ^ and goto El 

4. if Q nA (0,n-l)nA+(i)?) ^ 0, then picky,?/ G Q n>1 (0,n-1) such that d(y,y') = 1, 
W{yj= -1, and y' G A+(^), set % = i + 1, ^(y) = +1, = T^x) 
Vie A \ {y} and goto [3J 

5. if Q ni i(0,n-l)nA+(^) 7^ 0, then picky,?/ G Q n>1 (0,n-1) such that d(y,y') = 1, 
^f(y) = -1, and y' G A+ty?), set i = i + 1, = +1, $V) = +1, CO*) = 
Til>i_i{x) for any x G A \ {y, y'} and goto [3] 

6. set i — i + 1, y — (0, n), ^(y) = +1, ^(x) = T-^"_ 1 (x) for any x G A \ {y} and 
goto El 

7. set k n = i, * n = and exit. 

In the following lemma we summarize the main properties of the paths H n and *ff n 
defined above. In particular in items [2] and [3] we state upper bounds on the energy levels 
reached by those paths. We show that the addition of a slice of pluses to a rectangle of 
pluses can result in a net increment of the energy only if the length of the added slice 
does not exceed the critical length A (see (12. lip . ( 13. 5ft . and (13.6jl ). 

Lemma 3.1 Let a G S be such that o~(x) = +1 for any x G £^2,2(0), consider the path Q a 
defined by Then 

1. for any n = 3,...,L the configuration ifj n is such that tp n (x) = +1 for all x G 
Qn,n-i{0) , for any n = 3, . . . , L the configuration £ n is such that C, n (x) = +1 for all 
x G Q n ,n(0), in particular £ L = +1 and E L = {£ L }/ 
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2. for any n = 2, . . . , L we have 



E(ip n+1 ) - E(C) < (8 - 4/m) V and $ S n < E(C) + 10 - 6/t 



(3.5) 



5. /or any n = 3, . . . , L we have 



E(C) ~ E(4> n ) < (8 - 4/m) V and < E(^ n ) + 10 - 6h 



(3.6) 



4- we have 



- E(a) < T - 16(2 - fc) 



(3.7) 



where we recall V has been defined in Ii2.12\) . 

Proof of Lemma \3.1\ Item [1] is an immediate consequence of the algorithmic definition of 
Qcr. The proof of item [2] is similar to the proof of item[3l 

Item [3 Let * n := {ffi, . . . , . . . , with fc n > jfe > 1, such that = ^C-i 

for i = 2, . . . , k and ^j? = T 2 ip%; note that by construction ip™ = ip n , = £ n , and 
k n — k < n. By using (I2.24p . we get 



for i = 1 . . . , k — 1. If & n = A;, then (I3.6P follows immediately from (13.81) . In the case 
fc n > fc + 1, we shall prove that 



<S> { ^ +1 ,...^ L} <E(^) + 10-6h and E^D - E{^) < (8 - 4/m) V (3.9) 
The bounds (13.61) will then follow from (13.81) and (13.91) . 



We are then left with the proof of (13.91) . which can be achieved by discussing the 
following three cases. 

Case 1. There exist y,y' E Q n ,i(0,n — 1) such that ^(y) = — 1, y' E A+(^£). The 
configuration is defined at the step 0] of the algorithm; it is immediate to see that 
all the configurations ip™, with i = k + 1, . . . , k n , are indeed defined at the step HI Then, 
by using (12.181) . see also figure [TJ we get the following bounds on the transition energies: 



< Em + 2(1 - h) and > ^(C+i) + 2(1 + h) (3.10) 

for any % = k, . . . , k n - 1. By using (I3~T0|) . (12351) . and we get 

$ W ^ +1 ,..,^} < + 2(1 - h) and £(^J - < ~±h{k n - k) (3.11) 



which, recalling k n > k + 1, imply (13.91) . 

Case 2. There exist G Q n ,i(0,n — 1) such that i/j%(y) = — 1, ^iv') = +1> anc ^ 
^(^fc) ^ Qn,i(0, n — 1) = 0. The configuration t/^+i is defined at the step [5] of the 
algorithm; it is immediate to remark that all the configurations iftf, with % — k + 1, . . . , k n , 
are instead defined at the step HI 



%,«= £ W) and £(C)>£(C+i) 



(3.8) 
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Let i = k, . . . ,k n — 1, let y, y' G Q n ,i(0, n — 1) be the two sites which are picked up 
by the algorithm, let A, be the collection of the sites in Q n ,i(0, n — 1) different from y, y' 
and such that they become stable plus sites at this step of the path; more precisely, let 
Aj := A^^JYi) \ (^(V'?) y'})- Note that the update of the sites in Aj has no energy 
cost since they follow the zero temperature dynamics T. 

By using (12.181) , see also figure [fl we get the estimates 

E(r k ,r k+1 )<Em+^i-h) E(r k+1 ,r k )>E(r k+1 ) + 2(i + h)(i + \A k \) 

<Em + 2(1 -h) >*W + i) + 2(1 + / l )(l + |A i |) 

(3.12) 

for any % — k + 1, . . . , k n — 1. If | Aj| = for any % — k, . . . , k n — 1, then it must necessarily 
be k n — k = n — 1. We get 

®W,i>i +1 ,.,nj <E(M)+4(l-h) and £(^J -E(ift) < 2-2h-Ah(n- 1) (3.13) 

Noted that 8 - 4/m = 2 - 2/i - Ah(n - 1) + (6 - 2h), the bound ([33]) follows since h < 3. 
Suppose, finally, that there exists i G {k, . . . , k n — 1} such that |Aj| ^ 0; hence 

^ { ^^ k+1 ,...^ k j<E(r k )+^l-h) and tftyJJ - < -±h(k n - k + 1) (3.14) 

Recall k n > k + 1; the bounds (13.141) imply (13.91) trivially. 

Case 3. For each y G Q n ,i(0, n — 1) we have ip k (y) — ~ 1- I n this case k n — k = n, is 
defined at the step El "0^+2 i s defined either at the step H] or at the stepEJ and ip% +i , with 
i — 3, . . . , k n , are defined at the step 0] of the algorithm. By using (12.181) . see also figure [fl 
we get 

E(M, ^ + i) < E(^) + 2(3 - h) E(^ +1 , > 

£« +1 ,V£ +2 ) < ^+i)+4(l-/ i ) E(M +2 ,M +1 ) > £(^ +2 )+2(l + fc) (3.15) 
E(W, VT+i) < + 2(1 - &) > E(^ +1 ) + 2(1 + h) 

for % — k + 2, . . . , k n — 1; see the Figure for a graphical representation of the estimates 
(I3.15p . Note that the equalities hold, for instance, in the case ip^(x) = —1 for any 
x G 0Q„,i(O,n - 1) \ Qn,„-i(0). By using (13151) . fT2T26|) . and (I2TT9D we get 

<& {VW+1 ,...,^ } < + 10-6/1 and £(Vt) - (V£) < [8 - 4/i(fc n — k)]=S — Ahn 

(3.16) 

which imply (13. 9p . 

We remark that in this case 3 the path {ip k , ■ ■ ■ , ip kn } realizes the standard growth of 
the rectangular plus droplet ip n up to the square droplet ^ via the formation of a unit 
plus protuberance in the slice adjacent to one of the longer sides of the rectangle and the 
bootstrap percolation plus filling of the same slice. 

Item|H Let 77,77' two consecutive configurations of the path Q a , we shall prove that 

E(r], rf) - E(a) <Y — 16(2 - h) (3.17) 

The bound (13.71) will then follow, see (I2.26p . Recall the critical length A has been defined 
in (I2.1ip and consider the following four cases. 
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2(1 + h) 




Figure 5: Graphical representation of the estimates (I3.15p . 

Case 1. The configurations 77,77' belong to H n for some n < A — 1. This case is similar to 
the case 2. 

Case 2. The configurations 77,77' belong to ^ n for some n < A. By using (12.261) . (13. 5p . and 
(13. 6p we have 



£(77,7/) < < #(> n ) + 10-6/i 

< Eu> n ) - E^ 1 ) + w 3 ) + w 3 ) - ^(e 2 ) + ^(e 2 ) 

+ 10-6/i 

n-l A-l 

< £(cr) + 18-14/i + 8 J^[2 - < E(o-) + 18 - Uh + 8 J^[2 - /ii] 

j=3 i=3 

where we have used that 2 — hi > for z < A — 1 and £ 2 = a. The bound (13.171) follows 
easily. 

Case 3. The configurations 77, 77' belong to S n for some 77 > A. Note that for n > A the 
bounds (13.51) and (13.61) on the differences of energy become trivial since 8 — 4/777, < 0, hence 
< E(ifj x ). Then 

E(r), r/') < $ H n < E(C) + 10 - 6/7 < E(^ x ) + IO-6/7 

where in the first inequality we used (I2.26p . in the second the bound (I3.5p . and in the last 
the fact that E(£ n ) < E(t/j x ). To get (fXTTj) we then perform the same computation as in 
the case 2. 

Case 4. The configurations 77,77' belong to ty n for some n > A + 1. This case is similar to 
the case 3. □ 

Proof of item[J\ of Theorem \2.3[ Let a G S \ {—I}. If er = +1 the statement of the lemma 
is trivial; we then suppose a 7^ +1. Since by hypothesis a 7^ —1, there exists x G A such 
that o~{x) = +1; without loss of generality we suppose cr(0) = +1. Consider the path 
uj := {a, a 1 , a 2 , a 3 } with 
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- cr 1 is such that <r 1 (x) = +1 for all x G <52,i(0) and a l (x) = Ta(x) for all x G 

a\Q 2 ,i(o); 

- a 2 is such that a 2 (x) = +1 for all x G <5 2 ,i(0) U <5i(0, 1) and o" 2 (x) = Ta 1 {x) for all 
i£A\[Q 2)1 (0)UQ 1 (0 ) l)]; 

- cr 3 is such that a 3 (x) = +1 for all x G (52,2(0) and <7 3 (x) = Ta 2 (x) for all x G 

A\g 2 , 2 (o). 

By definition the path uo + J1 CT 3 starts at a and ends in +1_, i.e., u + fi CT 3 G 0(<7, +1), 
moreover we shall prove that 

< + T (3.18) 

The item [1] of Theorem 12.31 will then follow. 

To prove (13.181) we first consider the path uj\ by using (12.181) . see also figure [U we get 
the following bounds on the transition energies: 

E^a 1 ) < E(a) + 2 • 2(3 - h) E(a\a) > E^a 1 ) 

E(a\a 2 ) < Eia 1 ) + 2-2(1 - h) + 2(3- h) E(a 2 , a 1 ) > E(a 2 ) (3.19) 

E(a 2 , a 3 ) < E(a 2 ) + 3 ■ 2(1 - h) E(a 3 , a 2 ) > E(a 3 ) + 2(1 + h) 

By using fl3+9|) . fl2T26l) . fl2+9l) . fl2+2l . and the definition fl23TD of the critical length A, it 
is easy to show that 

$u, - E(a) < 28 - 16h< T (3.20) 

and 

E(a 3 ) - E{a) < 26 - ISh (3.21) 
We consider, now, the path S1 CT 3 ; by using (13.71) and (13. 21 j) . we get 

$n CT 3-£(<x) = $n a3 -E(a 3 )+E(a 3 )-E(a) < r-16(2-/i)+26-18/i = T-2{3+h) (3.22) 

The inequality fl3~T8l) follows from (13~20|) and (l3~22l) . □ 

5.^. The variational problem 

Item [2] of Theorem 12 . 31 deals with the determination of the minimal energy barrier between 
the metastable state —1 and the stable one +1, more precisely with the computation 
of $(— 1, +1). In the context of serial Glauber dynamics this problem has been faced 
with different approaches each suited to the model under exam, see [OV] and [MNOS, 
Section 4.2], where a quite general technique is described. All these methods rely on the 
continuity of the dynamics, namely, on the property that at each step only one spin is 
updated. 

In the case of parallel dynamics, see [CN], the lacking of continuity increases the 
difficulty of the computation of the communication energy between the metastable and 
the stable state. We follow, here, the method proposed in [CN] which is based on the 
construction of a set Q C 5 containing —1, but not +1, and on the evaluation of the 
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transition energy for all the possible transitions from the interior to the exterior of such 
a set Q. 

To define the set Q we need to introduce the two mappings A, B : S — > S. Let a G S, 
we set Aa := a if E(a x ) > E(a) for any x G A+(er), otherwise Ao := a x where x is the 
first element of A+ (a) in lexicographic order. The map A flips the first, in lexicographic 
order, unstable plus spin of a to which corresponds a decrease of the energy. Under the 
effect of the map A the number of pluses decreases, but only unstable pluses are flipped. 
Let a G S, the configuration Ba G S is such that for each x G A 

Ba(x):= {7^ if * eA i~» (3.23) 
I a[x) otherwise 

Note that the operator B performs a single step of bootstrap percolation; relatively to a, 
it flips all the minus unstable spins and, among the stable minus spins, only those with 
two neighboring minuses. 

In the sequel a relevant role will be played by the configuration B Aa, for any a G S; 
recall the definition of fixed point of a map given at the end of Section 12. 21 The sole 
unstable positive spins in Aa are those corresponding to energy increasing flips. Starting 
from Aa, the map B, which flips the minus spins with at least two plus spins among 
the nearest neighbors, is applied iteratively until a fixed point is reached. It is easy to 
show that the pluses in such a fixed point form well separated rectangles or stripes winding 
around the torus; more precisely, the pluses in B Aa occupy the region |J™ =1 Qe^^^i) c 
A, where n, £\^, £i j2 , • • • , £ n ,i, £n,2 are positive integers and Xi G A for any i = 1, . . . ,n, with 
Qi il ,£ i2 (xi) being pairwise not interacting (see Section [2TTT) . Note that, depending on the 
values of £ it i, £ it2 , the set u £ i2 (xi) can be either a rectangle or a stripe winding around 
the torus. 

We can now define the set Q. Let a G S, consider B Aa, and, provided B Aa ^ — 1, 
denote by Qe ilt e i2 {xi) the collection of pairwise not interacting rectangles (or stripes) 
obtained by collecting all the sites y G A such that B Aa(y) = +1. We say that a G Q 
if and only if B Aa = — 1 or minj^ x, < A — 1 and maxjfq,^^} < L — 2 for any 
i = 1, . . . ,n. Note that configurations a such that B Aa contains plus stripes winding 
around the torus A do not belong to Q. 

In general Ta ^ B Aa, this means that B Aa is not necessarily the result of the zero 
temperature dynamics started at a. This is not a problem when looking for the minimal 
energy barrier between — 1 and +1, provided the energy of such configurations is larger 
than T. The definition of Q is indeed satisfactory because we can prove the following 
Proposition 13.21 on which the proof of items [2] and [3] of Theorem 12.31 is mostly based. To 
state the lemma we need one more definition: recall the set C is defined as the collection 
of configurations with all the spins equal to —1 excepted those in a rectangle of sides A — 1 
and A and in a pair of neighboring sites adjacent to one of the longer sides of the rectangle. 
Then, given 7 G C, we let 7r(7) C S the set whose elements are the two configurations that 
can be obtained from 7 by flipping one of the two plus spins in the pair attached to one 
of the longer sides of the plus spin A x (A — 1) rectangle. We also let V be the collection 
of all the configurations with all the spins equal to —1 excepted those in a rectangle of 
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sides A — 1 and A and in a single site adjacent to one of the longer sides of the rectangle. 
Finally, we let TZ be the collection of rectangular droplets with sides A — 1 and A. By 
using ( 12.251) . we have 

E(K) - E{-1) = -Ah\ 2 + 4/iA + 16A - 8 = T - 10 + 6h (3.24) 

where we have used in the last equality the definition (12 . 1 2 j) of T. By using f)2.13p . we 
have the easy bound 

E(TZ) - E(-l) < 8A + 4J% (3.25) 

Proposition 3.2 With the definitions above, for h > small enough and L = L(h) large 
enough, we have 

i- -I e g, +i e s \ g, and c c s \ g ; 

2. for each r\ £ g and ( £ S \ g we have E(r), () > E(—V) + T; 

3. for each t] £ g and £ £ S \Q we have E(rj, Q = E(—l_) + T i/ and only if Q E C and 
V e AO- 
Proof of item\M of Theorem \2.3\ Since —1 £ g and +1 £ S\g, see item[T]in Proposition ^. 2\ 
we have that any path a> = {uj\, . . . , u n } such that U\ = —1 and u n = +1 must necessarily 
contain a transition from Q to S \ g, i.e., there must be i £ {2, . . . , n} such that Ui-i £ g 
and uii £ S \ g. Thus, item |5] in Proposition 13.21 implies that $ w > E(—l) + T; since the 
path uj is arbitrary, it follows that 

H-L+l)>E(-l) + r (3.26) 

To complete the proof of (I2.29P we need to exhibit a path connecting — 1 to +1 such 
that the height along such a path is less than or equal to E(—l) + T. Consider the path 
uj := { — 1, a 1 , a 2 , a 3 , a 4 } with a 1 the configuration with all the spins equal to minus one 
excepted the one at the origin, a 2 the configuration with all the spins equal to minus one 
excepted the ones associated to the sites in the rectangle Q2,i(0), o~ 3 the configuration 
with all the spins equal to minus one excepted the ones associated to the sites in Q2,i(0) U 
Qi(0, 1), and a 4 the configuration with all the spins equal to minus one excepted the ones 
associated to the sites in the square (^2(0). 

By definition, the path u) + Q a 4. starts at —1 and ends in +1, i.e., u + Q a 4 £ 0(— 1, +1). 
Moreover, we shall prove that 

$ w+ n CT4 - E(-l) < T (3.27) 

The inequality fl3T27l) . together with (I3T26D . implies fl2T29l) . 

We are then left with the proof of (I3.27p . We first consider the path uj; by using (12.18!) . 
see also figure [U we get 

E(-l, a 1 ) = E(-l) + 2(5 - h) E(a\ -1) = E(a l ) 

E(a\a 2 ) =E(a 1 ) + 2 ■ 2(3 - h) E{a 2 ,^) = E{a 2 ) + 2(1 - h) 

E(a 2 ,a 3 ) = E{a 2 ) + 2-2(1 - h) + 2(3- h) E{a 3 , a 2 ) = E{a 3 ) + 2(1 - h) 1 > 

Eia 3 ,^) = E(a 3 ) + 3 ■ 2(1 - h) E{a\a 3 ) = E{^) + 2(1 + h) 
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-1 



Figure 6: Energy landscape for the path oj. 



see figure M for a graphical representation. 

By using (jg^gjl . (G^5j) . fICT]) . and the definition (121TD of the critical length 

A, it is easy to show that, provided h is chosen smaller than 3 + y/E, 

$ w - E{-1) < 34 - Uh< T (3.29) 

and 

E(a A ) - E{-1) < 32 - 16h (3.30) 
We consider, now, the path S1 CT 4 ; by using (13.71) and (I3.30p . we get 

- E(-l) = $ n<j4 - E(a A ) + £(a 4 ) - E(-l) < T - 16(2 - h) + 32 - 16/i = T (3.31) 

The inequality (13.271) follows from ( 13.291) and ( 13.311) . This completes the proof of item [2] 
of Theorem 12.31 □ 

Proof of item\^ of Theorem \2.3[ The item follows from item [2] of Theorem 12.31 and item [3] 
of Proposition 13.21 □ 



4. Proof of Proposition 13.21 

In Section S3] we shall prove Proposition 13.21 concerning the solution of the minmax prob- 
lem. Some preliminary lemmata are stated in advance. More precisely, in Section 14.11 
we state the Lemmata I4.1H4.4I concerning energy estimates for the maps A and B (see 
Section 13.21) . In Section H~2l the Lemmata 14.51 and 14.61 concerning properties of rectan- 
gular droplets (see Section |2~5|) . are stated. Section I4T31 is devoted to the comparison of 
configurations in Q (see Section I3T21) and in Q c . 
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4-1. Energy estimates for the maps A and B 

In Lemma I4~T1 we give estimates on the energy of the configurations obtained by applying 
the maps A and B. For any a £ S we let 

Na(o-) := ]T[1 - 5 a{x) ^ {x) ] and N B (a) : = ^[1 - &Aa( x ),BA*( x )} (4-1) 

with 5 the Kronecker 5. Note that Na{o~) is the number of plus spins which are flipped 
by the iterative application of the map A to a, while Nb(<j) is the number of minus spins 
which are flipped by the iterative application of the bootstrap percolation map B to Aa. 

Lemma 4.1 Let a £ S and h > small enough. Then 

1. we have 

E(a) > E(Aa) + (2 - 10h)N A (a) (4.2) 

2. we have 

E(Ao) > E(BAcr) + AhN B (a) (4.3) 

In order to prove Lemma 14.11 we state Lemma 14.21 on some properties of unstable 
plus spins and Lemma 14.31 concerning an energy estimate for a single application of the 
bootstrap percolation map B. Recall (13. 1H . (13. 21) . and ( 13.31) : recall also that, given a £ S 
and x £ A, the configuration o~ x has been defined in Section 12.21 as the one obtained by 
flipping the spin of a associated with the site x. 

Lemma 4.2 Let a £ S; for h > small enough, we have that the following statements 
hold true: 

1. if there exists x £ A+(<r) such that E(o~ x ) > E(a), then \d{x} fl A~(c)| < 1, i.e., 
there exists at most one nearest neighbor of x which is stable w.r.t. a and such that 
the associated spin is minus one; 

2. if there exists x £ A+(er) such that E(a x ) < E(a), then E(a) > E(a x ) + 2 - lOh; 

3. if E(a x ) > E{o) for any x £ A+(cr), then 

2\KUo-)\ + 3|Ais(ff)| < Z\K x {<t)\ + 4|A; 3 (a)| (4.4) 

Proof of Lemma \4-2\ Let x £ A+(<r), then a(x) = +1, a x (x) = —1, and S a (x) < 0; by 
using ( 12.91) . we get 

E(a x ) - E(cr) = 2h-2+ ^ {\S a {y) + h\ - \S a (y) - 2 + h\) (4.5) 

y£d{x} 

Note that, since a(x) = +1, we have that S a (y), with y £ <9{x}, can assume the values 
—3, —1, +1, +3, +5; by performing the direct computations one shows that 

\S a (y) + h\- \S a (y) -2 + h\E {-2, 2h, +2} (4.6) 
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for y G d{x}. 

Item [fl Let x G A+(cr) such that ^(cr 21 ) > E'(er); since S CT (?/) < for y G <9{x} H A~(cr), 
by using ( 14.51) we get 

E{a x ) - £(<r) = 2h- 2(1 + |«9{x} n A b 7(a)|) + ^ fl&fo) + h\ - \S a (y) - 2 + h\) 

yed{x}\A-(a) 

Suppose, by the way of contradiction, that \d{x} fl A~(c)| > 2, then we have 
E(o*)-E{a)<2h-6 + ^ (\S a (y) + h\ - \S a (y) - 2 + h\) 

yed{x}\A~(a) 

By (14.61) we obtain \S a (y) + h\ — \S a (y) — 2 + h\ < 2 for y G <9{x}, and, noting that 
|<9{x} \ A" (a) | < 2, we finally get E(a x ) - E(a) <2h-6 + A = 2h-2<0, which is in 
contradiction with the hypothesis. 

Item El Let x G A+(cr) such that E(a x ) < E(a). Recalling (ED)]) and (j4~T)l) . we have 
that the number of sites y G <9{a;} such that \S a {y) + h\ — \S a {y) —2 + h\ = +2 must be at 
most equal to the number of sites y G d{x} such that \S a (y) + h\ — \S a (y) — 2 + h\ = —2, 
otherwise it would be E(a x ) — E(a) > 0. Thus, un upper bound to the sum in (14.51) is 
found when all the y G d{x} are such that \S a (y) + h\ — \S a (y) — 2 + h\ = 2h. We then 
get 

(\SM + h\- \S a (y) -2 + h\)< 2h\d{x}\ = 8h 

yed{x} 

from which E(a x ) - E{a) < -2 + Wh follows. 

Item[3l Consider o G S such that E(a x ) > E(o~) for any x G A+(cr) and let r a (y) = 1 
if y G A~ and r a (y) = otherwise. Recall that Al 5 (a) = and recall (13. 3ft ; by exploiting 
the first part of this lemma we get 

E E r <^) = E E r <^) + E E r <^) ^ 2 i A ±i( ff )i + 3 i A -3(^)i 

x£A+(a)y£d{x} xeAt^a) yed{x} x&A± 3 (a) yed{x} 

On the other hand, a site in A^ 1 (<r) is nearest neighbor of at most three sites in A+, 
indeed the number of unstable pluses neighboring such a site can be less than three since 
some of the pluses can be stable ones, and a site in A^ 3 (<t) is nearest neighbor of at most 
four sites in A+; then we have 

E E ^(y) < 3\A~ +1 (a)\ + 4\A^(a)\ 

x£A+(a) V&d{x} 

The inequality (14.41) follows trivially from the two bounds above. □ 

Lemma 4.3 Suppose h > small enough. Let a G S, suppose E(a x ) > E(a) for any 
x G A^(cr). Then 

E(a) > E(Ba) + 4/i|A>_ 1 (ct)| (4.7) 
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Recall that A > _ 1 (a) is exactly the set of sites whose associated spin flips under the action 
of the bootstrap percolation map B (see (13.231) ). 

Proof of Lemma \4.3\ To compare E(o~) and E(Ba) we shall use (I2.19p and suitable bounds 
on E(a, Ba) and E(Ba, a). Recall (I2.18p . see also figure HJ and the definition (I3.23|) of the 
bootstrap percolation map B; we have that in the forward transition from a to Ba the 
energy costs are those associated to the flip of the stable minuses with two neighboring 
pluses and those associated to the permanence of the unstable pluses. More precisely, we 
have 

E(a, Ba) = E(a) + 2(1 - h^AZ^a)] + 2(1 - h^AZ^a)] + 2(3 - h)\A± 3 (a)\ (4.8) 

On the other hand in the backward transition from Ba to a the energy costs that must 
be surely paid are those associated to the reverse flipping of the pluses that have been 
created in the forward transition; more precisely, we have 

E(Ba, a) > E(Ba) + 2(1 + h)\AZi(<r)\ + 2(3 + h)\A^{<T)\ + 2(5 + h)\AZ 3 (a)\ (4.9) 

Note that in (14.91) it is not possible to take advantage from the permanence of the possible 
unstable pluses in Ba, because, as we shall see in the proof of item [2] of the Lemma [4.11 
we have A+ (Ba) = 0. 

To complete the proof we have to distinguish two cases. Suppose, first, that AZi(a) = 
AZ 3 (a) = 0; by using flM}, flU), and flCTD . we get 

E(a) > E(Ba) + 4/i|A: 1 ( ( j)| - 2(3 - h)\A± 3 (a)\ + 2(3 + h^Az^a)] 

The bound (14.71) follows noting that, in this case, A>_ 1 (cr) = Al 1 (cr) U A^ 1 (cr) and (14.41) 
reduces to |A^ 3 (er)| < |A^ 1 (a)|. Suppose, now, that either AZ x (a) ^ or AZ 3 (a) ^ 0. By 
using (14.41) we have 

3|Al,( ff )| < 2|A± 1 (^)| + 3|A+,(a)| 

< 31^(^)1 +MA- 3 (a)\ < 3|A; x (a)| +5\AZ 3 (a)\ 

Since either |Al 1 (a)| > or |A7_ 3 (cr)| > 0, we have that 

\AtMl + 3|A+s(a)| < SIA;^)! + 5|A; 3 (a)| (4.10) 

We shall prove that, provided h < 1, 

(1 - /OlAi^)! + (3 - ^)|A+3(tr)| < (3 - h^A-^l + (5 - h)\AZ 3 {a)\ (4.11) 

First of all we note that the inequality (14.111) is equivalent to 

h[\AZ 3 (a)\ + ^(^1-^(^7)1-^3^)1] 

< |A; 3 (a)| + lA;,^)! - |A+!(a)| - \A± 3 (a)\ + [2^(01 + 4|A; 3 (a)| - 2|A+ S (a)|] 

which is trivially satisfied when the left hand side is negative or equal to zero, since (14.101) 
implies that the right hand side is strictly positive; on the other hand, if the left hand 
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side is strictly positive, recalling that h < 1, the inequality will follow once we shall have 
proved that 2|A+ 1 (<t)| + 4|A+ 3 («r)| - 2\A%(a)\ > 0. 

To get this last bound we note that by using item [T] in Lemma 14.21 it follows that for 
each site belonging to Al 3 , there are at least three unstable minus spins among the four 
nearest neighboring ones. Hence, we get |A~(<r)| > (4/3)|Al 3 (cr)|. Moreover, noted that 
| A- (a) | = | A^ (a) | + |A; 3 ((t)|, we also get 

2\A- +1 (a)\ + A\A- 3 (a)\ - 2\At 3 (a)\ 

= 2|A; 3 (a)| +2[|A- 1 ( C t)| + \A~ 3 (a)\ - \A± 3 (a)\) 

= 2|A;>)| + 2[|A-(^)|-|Al 3 (a)|] >2[|A; 3 ( C 7)| + l/3|Al 3 (a)|] >0 

Finally, the bound (14 .TP follows easily by using ( 12. 19ft . (14.81) . ( 14.91) . and the inequality 
@H). □ 



Proof of Lemma \4-l\ Item [TJ The bound (14.21) is proven easily by applying iteratively 



item [2] of Lemma 14.21 

Item [21 Suppose B Aa = B n Aa for some integer n. We first note that by Lemma [4.21 
each site x G A^(Aa) has at least two neighboring minuses which are unstable w.r.t. 
Aa, more precisely \d{x} fl A~(Ar)| > 2. Recall the definition (13.231) of the bootstrap 
percolation map B; since A~(Aa) C A>_ 1 (Acr), the minuses in d{x} fl A~(Aa) flip under 
the action of B. Hence, \d{x} n A + (B~Aa)\ > 2. We then have A+(BAa) = 0; in other 
words all the unstable pluses in Aa become stable after the application of a single step of 
the bootstrap percolation. 

By definition of the bootstrap percolation map we also have that A^(B l Aa) = for 
any i = 2, ...,n, i.e., no site in A + (B l Aa) is unstable w.r.t. B l Aa. Note, finally, that 
E((Aa) x ) > E(Aa) for any x G A^ (Aa). The theorem then follows by applying iteratively 
Lemma 14.31 □ 

Let cr g S, we refine the estimate (14. 2 p by considering the plus spins that are flipped 
by the iterative application of the map A and are associated with sites outside the support 
of the configuration B Aa. Let the branch of a be 

L(a) := \A+{a) \ A + (BAa)\ (4.12) 

i.e., the number of pluses outside the rectangles of B Aa which are flipped by the map A; 
note that L(a) < N A (a) (see ( H~Tj) ). 

Lemma 4.4 For any a G S such that L(a) > 1, we have that 

\ ,-./-r \ ^ f 6 — 2h if L(a) = 1 

E ^ ~ ( ^ > { l0 _ 6h+ (2 _ m{L{a) _ 2) lf L \J) > 2 ^ 4 - 13 ) 



Proof of Lemma\4-4\ Let a G S such that £(cr) = 1, the set A + (a)\A + (B Aa) has a unique 



element x. There exists a natural number j such that A 7 1 cx(x) = +1 and A^a(x) = —1. 
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For y G d{x} n_(A+(5 Ar)) c we have |5Ui ff (y) + /i| - iS^-i^y) + /i| = 2, while for 
y G <9{x} fl A + (£> At) we have the trivial bound \S^j a (y) + h\ — \S^j-i a (y) + h\ > —2. 
Since \d{x} n (A + (5lff)) c | > 3 and \d{x} n A+(5 Ar)| < 1, by using (Q we get 

E(A^a)-E(A^a)=2-2h + £ (|«W</) + fc| - |S^-v(</) + 

yG9{x}n(A+(BA<j))c 

+ (\Sa>M +h\- \S Aj -i a {y) + h\) > 6 - 2h 

y£d{x}llA+(BAcr) 

(4.14) 

Recall, finally, that by definition the map A decreases the energy; then, by (14.141) . we 
have 

E(a) > E(A j - l a) > E(A j a) - 2h + 6 > E(Aa) - 2h + 6 

and the bound (14.131) follows. 

Let now a G S such that L(a) = 2; the set A + (er) \ A + (B Aa) has two elements x, y. 
Since B Aa = BAa y and L(a) = 2, we have L(a y ) = 1; by using Acr^ = Aa and (jQ5]l in 
the already proven case we have that 

E(a) - E(Aa) = E(a) - E(a y ) + E(a y ) - E(Aa) > E(a) - E(a y ) + 6 - 2h (4.15) 

In order to bound E(a) — E(a y ), we first note that by (12.91) we get 

E(a) - E(a y ) = -2h - ^ (1^^) + h\ - \S aV (z) + h\) (4.16) 

We distinguish, now, two cases. We first suppose that x g" {y}, i.e., the two sites 
x and y are not nearest neighbors. It is easy to prove that —(\S a (y) + h\ — \S a v(y) + 
h\) = +2. Moreover, note that the contribution to the sum (14.161) of all the sites in 
d{y}n(A + (B Aa)) c is equal to +2 excepted for at most one site whose contribution is equal 
to -2h. Note, also, that \d{y} n (A + (BAa)) c \ > 3; hence, we have that E(a) - E(a y ) > 
-2/i+(2-2/i) + 2 + 2-2/i-2, where the contribution of the site d{y}nA + (BAa), which 
possibly exists, has been bounded trivially by —2. The bound (I4.13P follows immediately. 

Suppose, now, that x G {y}, i.e., the two sites x and y are adjacent. The only not 
trivial case, see Figure [7J is the one in which both the sites x and y are at distance one 
from the set A + (BAa). Since the plus spins associated to x and y are flipped by the 
iterative application of the map A to a, the spin associated to at least one of the two 
sites in d{x, y} fl A + (B Aa) is equal to —1, see Figured Without loss of generality we let 
d{y}f]A + (BAa) = {y'} and a(y') = -1. It is easy to prove that -(\S a (y) + h\-\S a y(y) + 
h\) = +2, -(\S a (x) + h\- \S aV (x) + h\) > -2h, -{\S a (y') + h\- \S a y(y') + h\) > -2, and 
-(\S a (z) + h\ - \S a y(z) + h\) = 2 for each z G d{y} \ \x,y'}. Hence, by using (1436]) we 
get 

E(a) - E(a y ) > -2h + 2- 2h -2 + 2 + 2 = 4- Ah (4.17) 

The bound fT4TT3l follows by (14371) and (14351) . 

Let, finally, a G S such that L{a) > 3. Let % a suitable integer such that L(A t cr) = 2. 
The bound (I4.13P follows easily by using the Lemma [4.11 and (I4.13P applied to A l o. □ 
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Figure 7: The three cases studied in the proof of the Lemma l4~4l on the left the not trivial 
one. 



4-2. Energy estimates for rectangular droplets 

We first state and prove the following Lemma on some simple geometrical properties of 
rectangles on the lattice. 

Lemma 4.5 Let Qi umi , for i — 1, ... ,71, be pairwise disjoint rectangles with sides U, G 
N \ {0}, such that li < rtii for i = 1, . . . ,n, and semi-perimeter p := + rrii). 



1. We have 



1 n 

-/>^kmi (4.18) 



4 

i=i 

2. If there exists a positive integer k such that £ { < k—1 and nii < k for all i = 1, . . . , n, 
we have 



^^<-^--^m, (4.19) 

3. If n > 2 and U > 2 then 



2 r 2 

i=l i=l 



\p 2 >J2kmi + p (4.20) 

i=l 

Proof of Lemma \4-5\ Item [1] we have 

Y \ n 1 n 1 n n n 

i=l i=l i=l i=l i=l 

Item |2j we have 

n " 1 1 n 1 71 1 n 1 71 

^2 ^i m i = 2 ^2 2 limi ~ 2 ~ l ^ mi + n ^^ fc - 2 k + - - mj 

i=l i=l i=l i=l i=l i=l 
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which implies ( 14.191) . Item[3j note that 



i=l i i=l 

n 



= jjy^jh + m if - 4^1;™, + y^(7» + rrii) (I j + mjfj 

i i=l ij^j 

1 " 

= 4 ( ~ m ^ 2 + + m ^ + 

i=l i^j 

4 n 

- 4^^' + m i) = ^ 



i=i 

where in the second step we used n > 2 and in the last step Zj Am, > 2. The bound (14.201) 
follows. □ 

We introduce the notion of semi-perimeter of a multi-rectangular droplet. Let n > 1 
and £1, mi, ... , £ n , "^n integers such that 2 < £1, mi, . . . ,£ n , m n < L — 2, a E S a n- 
rectangular droplet with sides £\, mi, . . . , £ n , m n , we let 

n 

p(a):=J2& + m ^ ( 4 - 21 ) 

i=l 

be the semi-perimeter of the multi-rectangular droplet a. 

Lemma 4.6 Lei £, m two integers such that 2<£<m<L — 2 and <t G 5 rectangular 
droplet with sides £ and m. If £ < A — 1, we have 

E(a) - E(-l) >8£>0 (4.22) 

U £ < A — 1 and m > A + 1, we /jave 

£(0-) - E(TZ) > Ah{l - S h ) > (4.23) 

where we recall 1Z has been defined above Proposition ^.^ and 5h below $2.11\) . 

Moreover, for n > 1 integer, for any n-rectangular droplet rj e S with sides 1 < £% < 
rrii such that £{ < A — 1 and rrii < A for i = 1, . . . , n, we have that 



1 n 

E (rj) - E{-1) > (4 - 2h)p( V ) + -J2 m * ( 4 - 24 ) 



2 

i=l 



Proof of Lemma \4-6] Suppose £ < A — 1: by using (12. 25ft we have -E'(a') — £'(— 1) = 
-Ah£m + 8(£ + m) = (8 - 4/ii) + 8£; since £ < A - 1, the lemma follows. Suppose £ < A - 1 
and m > A + 1, by using fl2~23D we have E{a) - E(TZ) = Ah(m - A) [(A -£)- 5 h ], which 
implies (Q3l) . 
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We finally prove (14.241) . Recall that by hypothesis £{ < A — 1 and nii < A for any 
% = 1, . . . , n; by definition of multi-rectangular droplets and by using (14.191) with k = A, 
we have 

A 1 n 

|A + (r/)|<-p(r/)--^m, (4.25) 

i=l 

Now, by using (12.251) . (I4.2ip . (I4.19p . and the fact that the support of a multi-rectangular 
droplet is made of pairwise not interacting rectangles, we have that 

1 n 

E( V ) - E(-l) = -Ah\A + (r))\ + 8p(rj) > p(v)(8 - 2hX) + - J2 m * 

i=l 

which implies (14.240 since A < {2/h) + 1. □ 



4-3. Relations between configurations in Q and in Q c 

Consider a G Q and rj G Q c , in Lemma [4.71 we state a property relating the pluses in 77 to 
those in B Aa and we bound from below the transition rate A(<r, rj) (see (12. 201) ). 

Lemma 4.7 Let a G Q and t] ^ Q, 

1. we have 

\k + {rj) \ A + (BA(a))\ > 2 (4.26) 

2. we have 

A( CT>r ,)> ( ! 2 -f f f0r ^ (< r ) = ° (4.27) 

v ' /; - \ 4 - Ah for L(a) = 1 K J 

Proof of Lemma \4-7\ Item [TJ the item follows from the definition of the subcritical set 
Q. Indeed, if |A + (?7) \ A + (B A(a))\ < 1, we have that under the map A the positive spin 
outside A + (BAcr) is flipped, so that A + (B~Ar]) C A + (BAa). Hence 77 G Q, that is a 
contradiction. 

Item [2 from fEOHl) we get 

A^r?) = 2 ^ >2 ^ |S ff (z) + /i| (4.28) 

z6A: jEA\A+(Blir): 
iWW'X >7(z)S CT (2)<0 

If L(a) = 0, by KM . KM . the theorem follows. Indeed, in the r.h.s of KM there 
are at least two terms corresponding to sites x and y such that rj (x) = t](y) = 1, and 
S a (x) < —3,S a (y) < —3. If L(a) = 1, from (14.261) there exist two sites 

{x, y} C A+(r/) \ A + (BAa) (4.29) 

Note that, since L(a) = 1, one has S a (x) < —1 and S a (y) < —1. From (I4.28P we have 
the bound 

A(<7,77) > 2(1 - /i) + 2(1 - /i) (4.30) 
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and the theorem follows, see also Figure [H 



□ 



4-4- Proof of the Proposition \3.2\ 

Let a G S and suppose B Aa ^ — 1, there exist n{a) G N\ {0}, £i(a), rrii(a) integers larger 
than 2, and Xi{a) G A for i = 1, . . . , n(a) such that 



If B Aa = —1 we shall understand n(a) = 1, £i(cr) = mi (a) = 0, and p(a) = 0, see 
also (14.211) . Let a G S, we order the droplets in A + (B Aa) so that £i(a) A rrii(a) > A for 
i — 1, . . . , k(a) and £i(a) Arrii(a) < A — 1 for i = k(a) + 1, . . . , n(a); note that for a G Q we 
have k(a) = 0, while for a G Q c we have A;(cr) > 1. For the sake of simplicity, for a G Q° 
in the sequel we shall let r;(o") := ^,(cr) — A and qi(a) :— rrii(cr) — A for i — 1, . . . , &;(cr). 



Figure 8: Restricted sets on which we evaluate E(r], () in the proof of item [2] of Proposi- 
tion E2J 

Before starting the proof of the Proposition 13.21 we sketch the main idea. We shall 
define the subsets of the configuration space A$ C A4 C A3 C A% C A\ C Q, B2 C 
£>i C Q c , and reduce the proof to the computation of E(rj, () for 77 G A5 and (^82 (see 
Figure E>. We recall (1QT1L (P2|) . (Q]l . and let 

A := {o- G £ : £i(cr) V m^cr) < A for i = 1, ... , n(a)} A 4 : = {<r G A 3 : n(<r) = 1} 

A 2 := {a G Ai : p(o-) < 2A + 4, L(cr) < 4A + 42} A 5 := {cr G Aa : pO) = 2A - 1} 

A 3 := {a eA 2 : p(a) > 2A - 50} 



n(a) 



A + (B Aa) = |J Q k{(j)tmi{a) {xi{a)) 



i=l 




(4.31) 
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and 

B 1 := {a G Q c : £i(a), m^a) < L - 2 for % = 1, . . . , n(a)} 

k(r7) (A 321 

B 2 :={aeBx: AhN B (a)- Ah J2(ri(a) + qi (a) +r i (a)q i (a))< 10- 2h} 1 ' } 

i=i 

In order to bound E(rj, Q for 77 G Q and ( G Q c , we shall use the identity 

E( V X)-E(-1) = [E(rj)-E(A V )} + [E(A V )-E(BArj)} 

+[E(BArj)-E(-l)} + A(rj,0 

which is a straightforward consequence of the definition (I2.20p . Depending on the choice 
of rj, the different terms in the r.h.s. of the identity ( 14.331) will be properly bounded in 
order to get the theorem. 

Proof of Proposition \3.2i Item [TJ The proof is an immediate application of the definition 
of the set Q (see Section 13. 21) . 

Items [2j Step 1. Let t] G Q \ A\ and ( G Q c . There exists i G {1, . . . , n(r])} such that 
h(v) V mi(r)) > A + 1; hence, by using (Ol) . N B (rj) > and (jQgj) . we get 

£7(77, - £7(-l) > [£7(77) " E(Ar))} + " £(-!)] + Afo C) (4-34) 

where (14.221) and (I4.23P have been applied to each non-interacting droplet in B Aa to 
deduce that E(BAa) - E(K) > 0. Now, if L(r]) = 0, by using (Q3jl . (OL ^(77) > 0, 
P7| . and (13~24D . we get 

£7(77, C) - £7(-l) > [E{K) - £7(-l)] + 12 - 4/i > T 
On the other hand, if £(77) > 1, by using (f4T34j) . (f4TT3jl . (14^71) . and (|3T24j) . we get 

£7(77, C) - £7(-l) > 6 - 2h + [E(K) - £7(-l)] + 4 — 4/i > T 
Step 2. Let 77 G A! \ A 2 and C e £ c . By using fl433|) . and Nb(t]) > 0, we get 

£7(77, C) - £7(-l) > [£7(77) - E(A V )} + [E(BA V ) - £7(-l)] + Afa, C) (4.35) 

Now, suppose p(r)) > 2A + 5, by using fl435|) . A(?7,C) > 0, the definition fl2~TT|) . 

and (I2.13p . we get 

£7(77, C) - £7(-l) > (4 - 2h)(2X + 5) > 8A + 12 - Uh > T 

provided h > is chosen smaller than 1/6. Suppose, finally, L(rj) > 4A+43. If B Arj 7^ —1, 
by using (14.221) we get E(B Arj) — £7(— 1) > 0; note that this bound holds trivially also in 
the case B Aq = — 1. Hence, by using this bound, (14.351) . (I4.13p . and (12.131) . we get 

E(r], () - £7(-l) > 10 - 6h + (2 - 10/i)(4A + 43) > T 
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Step 3. Let 77 G A 2 and £ G £ c \#i. There exists i G {1, . . . , k(()} such that 4(C)Vmi(C) > 
L - 2. Since 77 G A we have that ^(77) < 2A + 4 and L(rj) < 4A + 42, then by using (14.191) 
with k = A we have |A + (r/)| < \A + (B Arf)\ +L(rj) < \p{rj)/2 + L(rj) < A 2 + 6A + 42. Given 
the magnetic field h > 0, the number of plus spins in 77 is bounded by a finite number; 
then we can choose L = L(h) so large that there exist an horizontal and a vertical stripe 
winding around the torus with arbitrarily large width and such that rj(x) is equal to — 1 
for each x in such two stripes. Since in B A(, there exists a rectangular droplet of pluses 
with one of the two side lengths larger or equal to L — 2; we choose L so large that 
A(rj,Q > r. By using, finally, (I2.20p we get E(rj,() — E( — l) > T, once we remark that 
E{rj) - E(-l) > E(BArj) - E(-l) > 0. 

Step 4- Let 77 G A? and ( G E\ \ £> 2 - By using Lemma |4~T1 and Na(C) > 0, we have the 
bound 

E(C) - E{-1) > E(BAQ - E{-1) + AhN B {Q (4.36) 
By (I2~25l) and (12321 it follows 

n(C) n(C) 

E(BA()-E(-1) = _4^53(A + r i (C))(A + ft(C)) + 85](2A + r,(C) + ft(C)) 

i=l i=l 
n(0 

= n(C)(r - 10 + 6/1) - ^(n(C) + ft(C))(4/iA - 8) 

8=1 

n(C) 

-4n(C)(/iA-2)-4/i5>(C)ft(C) 

i=i 

"(0 

> (r- 10 + 2/i) -4hJ2(ri(() + qi(0 +^(C)ft(0) (4.37) 

i=l 

where in the last inequality we used (12. lip and the fact that T > 10 — 6h. Hence, by 
rf436|) and (1437]) . we have 

"(0 

£(C) - £(-1) > T - (10 - 2h) + 4/iiV B (C) - 4/i 5>,(0 + ft(C) + ri(C)ft(C)) 

i=l 

Since C G #2 \ £1, we get £(C) - E(-l) > T. Finally, by the inequality in fl2TT9l . we get 
E(r),Q-E(-l)>r. 

Step 5. Let 77 G A 2 \ A3 and C G B 2 . We note now that £(5 Arj) - E(-l) > 0, which is 
trivial if B At] = —1, otherwise it follows immediately from (14.221) . By using this bound, 
fl4~33]) . Lemma SH N A (r]) > 0, and N B (r]) > 0, we get 

E( V ,0-E(-l)>A(riX) (4.38) 

We find, now, a lower bound to A (77, Q by multiplying the the minimum quantum 
2(1 — h), see Figured! times the number of flips against the drift in the transition from 
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77 to C- More precisely, 



A(t7, C) > 2(1 - h) \{x G A : rj(x)S v (x) > 0, r)(x)((x) < 0}| 

> 2(l-/0(|A+(C)|-|Afo)|) (4.39) 



with 



A(r)) := A + (B Arj) U A+(rj) \ A+(B Arj) (4.40) 

where we recall the closure of a subset of the lattice has been defined in Section 12. 1L 

Recalling that the application of the map A does not add pluses, the number of plus 
spins in the configuration ( can be bounded from below by adding the number of pluses 
in B A( to the branch L(Q of ( and subtracting the number of pluses N B (() added by 
the bootstrap map B. Namely, we have 

n(C) n(0 

|A + (C)| > n(C)\ 2 + A$>,(0 + ©(C)) + 5>(0ft(C) - ^(0 + HO 

1=1 1=1 
Now, by using that Q G B 2 , (I2.1ip . and L(() > 0, we get 



|A + (C)I > A 2 + ^(A(r,(C) + ©(0-^(0-^(0) -^^ + ^(0 



^ A ' - i A + 5> " ^'(O + ©(C)) > A 2 - |A (4.41) 
i=i 

where we also used A — 1 > 0. 

We next bound from above |A(t7)|. We first note that by using (I4.40j) and (I4.12p we 

get 

\Hv) | < \A + (BArj) \ + 5L(t?) (4.42) 

Now, suppose that A + (B Arj) ^ — 1; by using (14. 19}) with k = X and exploiting 77 G ^2? 
we conclude 

|A(t7) I < -\p{r,) + 20A + 210 (4.43) 

Suppose, on the other hand, that A + {B Aq) = — 1. By using H4.42j) . we get |A(?7)| < 
5L(r)) < 20A + 210; hence the bound (I4.43P holds since in this case p(rj) = 0. 

We finally bound A(rj,() by using the preliminary inequalities (I4.39p . (14.41 j) . and 
(033]); we have 

Afo, C) > 2(1 - h) [A 2 - ^A - l -Xp{rf) - 210] (4.44) 
Recall 77 G A 2 \ A 3 , then p(r]) < 2A - 51; hence by using (@~3ED, pi . and (12131) . we 

get 

1 53 

E( v ,0-E(-i)>r + -- Y + o(h)>r 

where in the last inequality we have chosen h > small enough. 
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Step 6. Let 77 G A z \ A± and C G #2- By using (Q3j| . LemmaHU ^(77) > 0, iV B (r^) > 0, 
driS]), and A (77, C) > 0, we get 

£(77, C) - E(-l) > E(BAr]) - E(-l) = -4h\A + (BAr])\ + p(rj) (4.45) 

Now, since 77 G A3 \ A4, we can use ( 14.201) to obtain 

E(r), C) - E(-l) > -h(p(r])) 2 + (Ah + 8)17(77) (4.46) 

Finally, by exploiting the properties of the parabola on the right-hand side of (I4.46P and 
recalling that for 77 G ^.3 \ A4 the semi-perimeter satisfies the bounds 2A — 50 < p(r]) < 
2A + 4, it is immediate to prove that the parabola attains its minimum at 77(77) = 2A — 50; 
hence, by using (14.461) and (12.121) . we get E(r], () — E(—l) > T for h > small enough. 

Step 7. Let rj G A4 \ A 5 and C e B 2 . By using (14331) . (I4T3D . and iV B (r/) > 0, we get the 
bound 

E( V , - E(-l) > [E(rj) - E(Ar])} + [E(B Ar,) - E(-l)] + A(r), () (4.47) 

Since rj G A4 \ A$, we have that n{rf) = 1 and then 2A — 50 < p{rf) < 2A — 2. We repeat, 
now, the same argument used at Step 6, but, since 71(77) = 1; we have to use (I4.18P instead 
of (OUj) : we then get 

E( V , C) - E(-l) > [E( V ) - E{Aq)] + [r - 10 + 0(h)] + A( V , () (4.48) 

Moreover, since 77.(77) = 1 and 79(77) < 2A — 2, by using the same arguments developed 
in the proof of (I4.26p . we get 

\A + (C)\A + (BArj)\ >3 (4.49) 

To complete the proof of the Step 7, we distinguish four cases by means of the param- 
eter L{rj). Consider, first, the case L{rj) > 3; by using (14.481) . (14.131) . and A(t7,£) > 0, it 
follows immediately E(rj, () — E(—l) > T. 

Consider, now, the case L(r]) = 2. We first note that by using (14.481) and (14.131) we 

get 

E(rj, - E(-l) > 10 - 6h + [r - 10 + 0(h)] + A(rj, () > T + A(r), () + 0(h) (4.50) 

The result £"(77, () — E(—l) > T will then be proven once we shall have obtained the 
bound A(t7,C) > 2(1 - h). 

To prove such a bound, we note that there exist x, y G A + (7/) \ A + (B At]); since 
x, y G A\A + (B Arj), by the definition of the two maps A and B, it follows that they cannot 
be both stable w.r.t. 77 (see Section[3]). If one of the two sites x and y, say x, is stable w.r.t. 
77, it is immediate to prove that x G dA + (B At]) and {y} = d{x} \ A + (B Arj). Since x 
and y are nearest neighbors, it follows that there exist no site in A~(r]) \ A + (B At]) which 
is unstable w.r.t. 77; hence, by using (I4.49p and (12.201) . it follows that A(r],Q > 2(1 — h). 
We consider, now, the case when both x and y are unstable w.r.t. 77. Suppose that either 
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((x) = +1 or C(y) = +1; from (l2~20l we have A(r],() > 2(1 - h). On the other hand, if 
((x) = ((y) = — 1, it is easy to see that, since L(rj) = 2, we have \A~(r]) \ A + (B Arj)\ < 1 
(recall the definition lEOjl ). Then, by using (jOSj) . it follows that A > 2(1 — h). 

Consider, now, the case L(rj) = 1. We first note that, by using ( 14. 48ft and (14. 131) . we 

get 

E( V , C) - E(-l) > 6 - 2h + [r - 10 + 0(h)] + A(rj, () > F - 4 + A(rj, C) + 0(/i) (4.51) 

The result E(rj,Q — E( — l) > F will then be proven once we shall have obtained the 
bound A(r/,C) > 3-2(l-/i). 

To prove such a bound we let a; the site such that {x} := A + (r/) \ A + (B At]). Suppose 
C(x) = +1; since x is unstable w.r.t. 77, by (I2.20p and (I4.49p . we have A (77, £) > 2(1 — h) + 
2(1 -/i) + 2(1 -h). On the other hand, suppose ((x) = -1; since |A~ (r/) \ A + (5 Ar/)| = 0, 
by it follows that A(t7, C) > 2(1 - h) + 2(1 - h) + 2(3 - /i). 

Consider, finally, the case L{rj) = 0. Recall (j2.20P ; the condition (I4.49P implies that 
A(?7, C) > 3 • 2(3 - h). Hence, by using also flPg|) . (jO]l . and ^(77) > 0, we get 

^(r/, C) - E(-l) > T - 10 + O(h) + 3 • 2(3 - h) > T 

Step 8. Let 77 e ^4 5 and ( G S2. We remark that, since p(r7) = 2 A — 1 and i\ V m\ < A, 
we have B ~M\ G ^. Hence, by using P~35]) . (jO]l . and (|5T53]> . we get 

^(r/, C) - > [^(r/) - E(A7])} + 4/i^(r/) + [r - 10 + 6/\] + A(t7, C) (4.52) 

To complete the proof of the Step 8, we distinguish four cases by means of the pa- 
rameter L{rj). Consider, first, the case L(rj) > 3; by using (14.521) . (14.131) . N B (rj) > 0, and 
A(?7, C) > 0, we get 

E{ri, C) - E(-l) > 12 - I6/7, + T - 10 + 6h > T 

Consider, now, the case L{rj) = 2; we let x, y be the two sites in A + (?7) \ A + (B Arj). 
We first note that, by using the inequalities (I4.52p and (I4.13p . we get the bound 

E( V , C) - E(-l) > 4hN B ( V ) + T + A(77, C) (4.53) 

Suppose, first, N B (r]) > 1; by (14.531) and A(rj,Q > 0, we immediately get E(r],Q — 
E(—l) > T. We are then left with the case N B (r)) = 0, i.e., A + (r]) D A + (B Arj); by using 
(I4.53p . the result E(r),Q — E(—l) > F will be proven once we shall have obtained the 
bound A(t7,C) > 2(1 - h). 

We note that |A+(?7) \ A + (_B ^.77) | < 1, indeed if by the way of contradiction x and y 
belonged both to A+(?7) \ A + (B Arj), then it should necessarily be x, y G dA + (B Arf) and 
d(x,y) = 1, namely, there would be a two-site protuberance added to the A x (A — 1) 
rectangle of pluses which is present in rj. Hence, we would have 77 G C C Q c , which is a 
contradiction. 

Suppose |A+(t7) \ A + (B Ai])\ = 1 and let x be the site in A+ (77) \ A + (B At]); since x 
is stable w.r.t. 77, we must necessarily have x G dA + (B Ar]) and y G d{x} \ A + (B Arj). 
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Note, also, that this implies | A~ (77) \ A + (B Ar])\ = 0. Thus, for ((y) = +1, in the sum 
in (I2.20p there is at least the term corresponding to y\ then we have A(?7, () > 2(1 — h). 
On the other hand, if it were ((y) = —1, recalling (14.261) we would have that in the sum 
in ( 12.201) there is at least a term corresponding to the flip of the spin associated with a 
site in A~(?7) which is stable w.r.t. 77, hence we would have A(t],Q > 2(1 — h). Suppose, 
finally, |A+(r/) \ A + (Blr?)| = 0; it is immediate to prove that [A~ (77) \ A + (Blr;)| < 1. 
Then we get A(t],() > 2(1 — h), since from (14.261) it follows that in the sum in (12.201) 
there is at least a term corresponding to the persistence of the spin associated with a site 
in A+(?7) or to the flip of the spin associated with a site in A - (77) which is stable w.r.t. rj. 

Consider, now, the case L{rj) = 1. We let x be the site in A + (?7) \A + (B At]), note that 
x in unstable w.r.t 77 and w is stable w.r.t. 77 for any w G A~(r]) \ A + (B Arj). We remark 
that by using (1452]) . (l4~T3l) . (IP) , and N B (r]) > 0, we get the bound 

E( V , C) - E(-l) > T - 4 + Ah + A( V , C) (4.54) 

and distinguish different cases depending on the number of plus spins in the configuration 
£ which are associated to sites outside the support of the configuration B At], that is on 
|A+(C) \ A+(BAr))\ > 2, see ( TC26T) . 

Suppose, first, |A + (£) \ A + (B Ar])\ > 3; since x G A+(7/) and w is stable w.r.t. 77 for 
any w G A~(r]) \ A + (B Arj), we have A(?7, () > 3 ■ 2(1 — h), for in the sum in (I2.20p there 
are at least three terms. 

We are left with the case |A+(C) \ A+(BAr)) \ = ^we let {y, z] := A+(C) \ A + {BAt]) 
and notice that it must be necessarily y,z G dA + (B Arj) and d(y, z) = 1, otherwise it 
would be ( G Q. Suppose, first, ((x) = —1; since x 7^ y, x 7^ z, and y and z are nearest 
neighbors, it follows that at most one of the two sites y and z is nearest neighbor of x. 
Then, since in the sum ( I2.20p there are at least two terms and one of them is greater or 
equal to 2(3 — h), we have A(?7, () > 2(1 — h) + 2(3 — h). By the previous inequality and 
(I4.54p we get E(r], £) — E{— 1) > V. Suppose, finally, ((x) = +1; without loss of generality 
we let y = x. Since z G d{x}ndA + (B Arj), by (12T201 we have A(r/,C) > 2(1 -h) +2(1 -h), 
with one of the two terms corresponding to the persistence of the plus spin associated to 
x in rj and the other corresponding to the flip of the minus spin associated to z in r\. By 
the previous inequality and (14.541) we get E(r], () — E( — l) > T. 

Consider, finally, the case L(r]) = 0. By using (T4T52D . g^D, N A (rf) > 0, N B (r]) > 0, 
and (HZ27D, we get 

E(r/, C) - E(-l) > T - 10 + 6h + 12 - Ah > T 

Item El Suppose ( G C and rj G vr(C), by using and frX20l) it follows E(rj, C) - 
E(-l) = r. 

Conversely, suppose 77 G Q and ( G Q c such that E(r],() — E(—l) = T. By using the 
results in the proof of the Item [2] above, see in particular the Step 8, we have that it must 
be necessarily 77 G A 5 , ( G B 2 , L(rj) = 1, |A + (C) \ A + (BArj)\ = 2, and ((x) = +1, with x 
such that {x} = A + (t]) \ A + (B At]), indeed for any different choice of rj and ( it has be 
proven E(rj, () — E(—l) > T. For configurations 77 and ( as above we have also proven 
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that BAr] ell, that A(r], C) > 2 • 2(1 - h), and that there exists z G d{x} n dA + (BAr]) 
such that ((z) = +1. 

Now, by using BAr] G TZ, A(t?, () > 2 • 2(1 - /i), fl433|) . fl4T3|) . (Ql) . and (13T24D . we 

get 

£(77, C) - E(-l) >6-2h + AhN B (n) + T - 10 + 6h + 2 ■ 2(1 - /i) = T + AhN B (rj) 

If it were Nb^) > 1 it would follow E(r],Q — E(—l) > T, then it must necessarily be 
N B (v) = 0- 

By the above characterization of 77 we have that rj G V\ then, by using (14.331) and the 
definition of the map A we get the following expression for the communication energy 
E(r),(): 

E(t), C) - E(-l) = 6 - 2h + T - 10 + 6h + A(r/, C) = T - 4 + 4/i + A (77, C) 

Since ((x) = ((z) = +1, we have that A(r/, () = 2 ■ 2(1 — /i) if and only if C(^) — +1 f° r 
all w e A + (BAr]). We then have that C G C and 77 G tt(C). □ 



A. Review of results in [MNOS] 

The proof of Theorem 12.11 is based on general results in [MNOS, Theorem 4.1, 4.9, and 
5.4] concerning the hitting time on the set of global minima of the energy for the chain 
started at a metastable state. We restate those results in our framework which is slightly 
different from the one considered in that paper (see the discussion at the beginning of 
Section £2]). 

Recall (12.271) . Let S s be the set of global minima of the energy (12.91) . For any a G S, 
let X a := {77 G S : E{j]) < E(a)} be the set of states with energy below E(a) and 
V a := 5>(<t, X a ) — E(a) be the stability level of a. Set V a := 00 if X a = 0. We define the set 
of metastable states S m := {r] G S : V v = max ffg5 \ S s V a }. We say that W(?7, C) C S is a 
gate for the transition from 77 G S to ( G S if and only if the two following conditions are 
satisfied: (1) for any a G W(?7, C) there exist a path to G 0(?7, (), such that $^ = $(77, C), 
and i G {2, . . . , |o;|} such that cjj = a and E(uii-i, cui) = $(77,^); (2) c<j PI W(t7, Q 7^ for 
any path u = 6(77, £) such that $ w = $(77, £). A function / : /3 G M — > f{f3) G M is called 
super- exponentially small (SES) in the limit (3 — > 00 if and only if limp^oo^l / j3) log f{(3) = 
—00. Given a G 5 and A C 5, finally, recall the definition of hitting time r\ given in 
(12. 4p and the notation E CT introduced just before it. 

Theorem A.l (restatement of Theorem 4.1 in [MNOS]) Let a G <S m ; for any § > 

0, there exist /3q > and K > such that, for any (3 > /3q, 

P ct (t£ s < e P v °-P & ) < e ~ K P (A.l) 

and 

P CT (r£ s > e ^ +/35 ) = SES (A.2) 



pca'croce.tex - 3 Luglio 2009 



39 



13:28 



Theorem A. 2 (restatement of Theorem 4.9 in [MNOS]) Let a G S m , then 

Jim - logE CT [r| 8 ] = V a (A.3) 

Theorem A.3 (restatement of Theorem 5.4 in [MNOS]) Let <j,r) G 5; consider a 
gate W /or £/ie transition from a to r\. Then there exists c > suc/i i/iai 

P.(r^ > t£) < e"^ (A.4) 

/or (3 large enough. 

The proof of Theorems IA. 1I - TA31 can be achieved by arguments much similar to the one 
developed in [MNOS]. To this purpose the main ingredients are the revised definition of 
cycle and the revised statement of [MNOS, Theorem 2.17] (see also [OV, Theorem 6.23]) 
and [MNOS, Theorem 3.1]. 

Let A C S, consider $(A, A c ), we say that A is a cycle if and only if 

max $(cr,r/) < $(A, A c ) (A. 5) 

cr,rj£A 

Let a G S; we say that the singleton {a} C 5 is a trivial cycle if and only if it is not a 
cycle. Given a cycle A C 5, we denote by -F(-A) the set of the minima of the energy in A, 
i.e., 

F(A) := {a E A: mm E(r]) = E (a)} (A.6) 

We also write E(F(A)) = E(a) with a G F(A). Noted that $(o, A c ) = $(«/, A c ) for any 
cr, a' G -F(A), we pick a G F(A) and set $(A) := $(o, A c ). 

Theorem A.4 (restatement of Theorem 2.17 [MNOS]) Let A C S be a cycle. For 
any a G A, r\ G A c , e, e' > 0, 5 G (0, e), and (3 > /arge enough 

P CT (rlc < e/W)-W))]+*. r^c = > e -/W*A))-*(A)]-*' (AJ) 

and 

P CT (rlc > efl*W-W))]-*) > 1 _ (A.8) 
Moreover, there exists k > stzc/i £/iat /or any cr, a' G A and /3 large enough 

P CT (r;, < r^ c ) > 1 - e"^. (A.9) 

Equation (1A.7I) is a bound from below to the probability that the chain exits a cycle A 
in a time smaller than the exponential of f3 times the height — E(F(A)) of the cycle 
plus e. In particular, it is stated that for such a probability the estimate is optimal when 
the exit from the cycle is achieved by touching a configuration rj such that $(77, A) is equal 
to $(A). In this case, for any e', provided /? is chosen large enough, such a probability 
is larger than exp{— j3e'}. Equation (I A. 81) is a bound from below to the probability that 
the chain exits a cycle A in a time larger than the exponential of (3 times the height 
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<&(A) — E(F(A)) of the cycle minus e; such a probability is larger than 1 — exp{— f3S} 
with 5 G (0,e). Equation (1A.9[) is a bound from below to the probability that the chain 
started at a G A visits another configuration a' belonging to the cycle A before exiting 
it; such a probability is larger than 1 — exp{— (3k,} with k > not depending on a and a'. 

Before stating the revised version of [MNOS, Theorem 3.1], we introduce the concept 
of metastable state at level V G R. We call metastable set at level V G R the set of all 
states with stability level strictly larger than V, i.e., Sy := {o~ E S : V a > V}. Any 
cr G Sy is such that for any path u starting from a and ending in a configuration with 
energy lower than E(a), the quantity $ w — E(a), that is the energy level reached along 
the path and measured with respect to a, is lower than V. 

Theorem A. 5 (restatement of Theorem 3.1 in [MNOS]) For any e > and (3 > 

large enough 

supP^rJ > e pv+pe ) = SES (A.10) 

Equation flA.lOh states that the probability that the chain started at o G S visits a 
configuration with metastability level V G R in a time larger than the exponential of (3 
times V plus e is super-exponentially small in f3. 

The proof of Theorem lA.5l can be achieved by repeating the same arguments developed 
in [MNOS] and based on [MNOS, Theorem 2.17]. The proof of Theorem IA.4I can be 
achieved by repeating the same arguments quoted in [MNOS] and developed in the proof 
of [OV, Theorem 6.23] (see also [OS, Proposition 3.7]). In particular, in the proof of (1A.8I) . 
a revised version of the so called reversibility property [OS, Lemma 3.1] is needed. 

Lemma A. 6 Let A C S be a cycle. For any a G F(A) and e > 0, there exist ft® > and 
c > 0, such that for any (3 > (3$, 

P CT (rlc < e ^)-E(F(A))]-^ < e -c P (A n) 

To prove Lemma IA.6I we develop an argument much similar to the one proposed by 
Olivieri and Vares to prove [OV, Lemma 6.22]. The difference between the two cases is 
in the fact that the Hamiltonian of the model ( 12. ip of the present paper depends on the 
inverse temperature f3, while [OV, Lemma 6.22] refers to a model (see [OV, Condition R 
in Chapter 6]) with Hamiltonian not depending on the temperature. Note, also, that the 
reversibility statement (1A.11I) is given in terms of energy-like quantities not depending 
on /3; on the other hand in the proof of the lemma the key property is the reversibility of 
the dynamics w.r.t. the Hamiltonian (12. 5p . 

Proof of Lemma \A.6l First of all we make explicit how hamitonian-like quantities differ 
from the corresponding energy-like quantities multiplied times (3. More precisely, given 
a, r] G S, we have 

H(a) = (3E(a)-^\og[l + e- 2 ^ {x)+h \] + |A| log 2 (A.12) 
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and 

H{a,Ti) = pE{<T,ri) + \A\ log2 (A.13) 

The equality flAl3l) follows easily from flAl2|) and (12T23D . We prove, now, (lA~T2l . Using 
( 12. 5p . (I2.9p . and recalling that coshx = cosh(— x) for any a; G R, we have that 

H(a)-(3E(a) = - ^ log cosh(/5|5 CT (x) + + (3 ^ + h\ 

- ^[logcosh^l^^) + h\) - logexp^l^^) + h\)\ 
xeA 

e f3\S<?(x)+h\ i e -j9|S ff (x)+/i| ]_ , e -2/3|5 CT (x)+fc| 

= ~ E log 2^l^>+^l = ~ E log 2 



which yields (1A.12[) 



Consider, now, an integer T > 2; recalling that the chain is reversible with respect to 
the Gibbs measure fi defined above (12. 5p . we have that 

T-l 

P CT (T^c < T) = fr>foO + E E P^Zl) ■ ■ ■p(tn-l,Zn)p{Zn,0] 

?eA c n=l £i,...,£ n eA 

= EW ff ^ + E E e-^»)-^(en,en-i)---Ke 1 ,^(en,e)] 

?eA c n=l f !,...,$„ 6 A 

(AAA) 

where the detailed balance ( I2.6P has been used to invert the order of the configurations in 
the terms p(a, £i), . . . ,p(£ n _i, £ n ). By using the definition (12.171) of transition Hamiltonian, 
we have that exp{-[i7(£ n ) - H (o)]} p(£ n , g) = exp{-[H(£ n , £) - H(a)]}. Noting, also, 
that 

Ci,-,fn-ieA 

we have 

P CT (r Ac < 20 = E + £ E e-^-O-^M P.K = £ n )] 

£eA c n=l § n eA 

£eA c n=l § n eA 

= e H M £ [e^^ + (T - 1) £ e"^] < T £ e"^ 
£eA c ceA CeA,feA c 



which yields 

P (7 (r yl c < T) < e H(a) T \A\ \A C \ max e- H{ ^ = e H{u) T \A\ \A C \ e -^^A,^ h(C,0 

CeA,£eA c 

(A.15) 



pca'croce.tex - 3 Luglio 2009 



12 



By using, finally, (TATL2]) . flAl3|) . and (EQ5l) . we get 



P.(t> < T) < eP E ^ T \A\ \A C \ e ~^ A ^ E «<$ exp { - ^ log [l + e -2/*|s„(*)+fc|] | 

(A.16) 

Since A is a cycle, we have that min^^g^c E((,l;) = <$>(A). Hence, remarked that 



< exp { - ^ log [1 + e-WW+M] } < (1 + exp{-2/?(5 + /i)}) HA| < (| 



3\|A| 



for /3 large enough, the bound (1A.8j) follows once we take T = e^*^ □ 
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